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O MÉTODO DOS GRADIANTES CONJUGADOS E A 
SUA APLICAÇÃO NA ANALISE DE ESTRUTURAS 


RESUMO 


Mostra-se que a utilização do método dos gradian- 
tes conjugados na resolução dos sistemas de equações 
associadas com a análise de estruturas torna possivel 
trabalhar directamente com as matrizes de rigidez dos 
elementos e, deste modo, evitar a formação da matriz 
de rigidez global. 


1. INTRODUÇÃO 


É sabido que, apesar de se dizer que a formulação 
matricial da Análise de Estruturas é especialmente ade- 
quada ao cálculo automático, essa formulação nunca é 
totalmente utilizada pelos programadores, 

Com efeito, em vez do chamado «stifíness method», 
estes utilizam o «direct stifness method», muito mais 
económico, em que a matriz de rigidez não se calcula 
por operações matriciais a partir das matrizes de rigidez 
dos elementos, mas por espalhamento [1] («scattering») 
& acumulação dessas matrizes sobre a matriz de rigidez 
global da estrutura. É portanto em termos desta técnica 
que se tem posto a problemáica da análise estrutural, 

Constitui objectivo do presente trabalho mostrar 
que outros processos podem ser utilizados porventura 
mais convenientes. 

O problema está intimamente relacionado com o da 
resolução de sistemas de equações a que a análise acaba 
por reduzir-se. Não admira portanto que o assunto da 
resolução de sistemas de equações seja abordado nas 
páginas seguintes. Tomar-se-ão como base os capítulos 
que Odent lhe dedica em [2]. 

Admite-se que o problema da resolução do sistema 
pode sempre converter-se no da minimização de um 
certo funcional. 
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E. R. DE ARANTES E OLIVEIRA 
Professor Catedrático do Instituto Superior Técnico 
L. F. GONÇALVES DE AZEVEDO 


Assistente do Instituto Superior Técnico 


SUMARY 


lt is shown that the use of the conjugated gradiant 
technique in the solution of the systems of equations 
associated with the analysis of struclures makes it 
possible to work directly with the construction of the 
global stiffness matrix. 


Lembra-se que, no caso dos sistemas de equações 
lineares da forma 
K q=0 (1.1) 


em que K é uma matriz simétrica e definida positiva, 
o funcional em questão é 


Flg='/2q'kKkq—Qa"a (1.2) 
Com efeito, diferenciando (1.2), obtém-se 
(1.3) 


5 F=ãqka—s qo = qiK q— O) 


im 


e, portanto, dado que dq é arbitrário, a estacionaridade 
de F é condição necessária a suficiente para a verifi- 
cação da equação. (1.1) 

Para mostrar que F é minimo no seu ponto de 
estacionaridade, há que determinar, diferenciando (1.3), 
a segunda variação de F. Obtêm-se 

Ss" F=5 qlKk q— 0) + ag!K sq (1.4) 
ou, em virtude de se tratar de um ponto de estaciona- 
ridade no qual (1.1) se verifica, 


8º F=a5q'K sq (1.5) 


Logo, visto que K é matriz definida prsitiva, 
 F> 0, o que revela tratar-se de um minimo. 

Admitir-se-á nas Secções seguintes que o sistema 
de equações é do tipo (1.1). Alguma coisa se dirá 
porém ao longo do texto sobre sistemas não lineares, 


2. MÉTODOS SEQUÊNCIAS 


Os métodos que vão referir-se são do tipo sequen- 
cial, isto é, tais que a solução ss determina a partir de 
uma sucessão de soluções aproximadas obtidas mediante 
uma rotina fixa de operações. 

A solução de ordem r+ 1, qrTi 
ordem r por uma fórmula do tipo 


é obtida da de 


a rti=q+N D' (2,1) 


em que o vector D' definee a direcção do incremento e 
NX é um escalar que fixa o seu valor. 

Na vizinhança de um ponto q o funcional F (g) 
é susceptível de desenvolvimento em série de Taylor, 
isto é, 


õ F 
la tm=F ta) + 2h; (5, + 
- - 7 Nog; fa 


1 5” F 
+ ob b ) (2.2) 
21;  !13q;êg;/a 


A fórmula (2.2) é exacta porque, sendo o funcio- 
nal quadrático, as suas derivadas de ordem superior à 
segunda anulam-se todas e os termos seguintes da série 
são portanto também todos nulos. 

Em virtude de (2.2) pode exprimir-se o valor de F 
no ponto q +' em função do seu valor e derivados no 


ponto q”. Atendendo a que o incremento é 


h; =N D, [2.3] 


vem 


à F 
Flg+H)=F (gq) +NI rr DE + 
- ô ú 


(NJ po 8? F 
ES ES D: D' (2.4) 


ou, designando por qg' (gradiante de F no ponto q) o 


ô F 
vector de componentes E ) e atendendo a que a 
dila 
5º F 
matriz de elementos é a matriz K, 
dq; dq; E 


7º 
FIO) =F (q) + QD + DK Dº 


ma O 


(2.5) 


Como determinar NX ? 

É evidente que convém N seja tal que o valor 
de F (gq'+') se aproxime tanto quanto possivel do mi- 
nimo, que corresponde à solução exacta. 

Nestas condições, não admira que o valor de À se- 
ja determinado derivando F (q'+!) em ordem a X e 
igualmente a zero a derivada, 


Ora, derivando ambos os membros de (2.5), 
obtém-se 
8 Flg'+>) 
= 9TD' + x DTK D'= O (2.6) 
A o E SME me 
o que dá 
gTD' 
fãs ua emma = (2.7) 
D'k D 


Como escolher agora os vectores D'? 

Várias técnicas podem ser utilizadas. 

Se os incrementos do vector das incógnitas se Ífize- 
rem sempre de modo que uma única incógnita seja alte- 
rada de cada vez, isto é, se todos os elementos dos vec- 
tores D' forem nulos excepto um, tem-se o chamado 


método univariante. 
So se fizer 


D' = —g (2.8) 
tem-se o chamado método da descida mais ingreme (*) 
(dado que o gradiante define a direcção da descida mais 
ingreme). 

Mostra-se que, no método da descida mais ingreme, 
a sucessão de soluções aproximadas converge para a 
solução exacta se a matriz K for definida positiva, como 


its 


se supõe, 

No entanto, no método da descida mais ingreme, 
a solução exacta é em geral atingida após um número 
infinito de passos. 

É portanto bastante mais vantajoso o método que 
vai passar a ser descrito, dito método dos gradiantes 
conjugados, no qual a solução exacta se atinge em n 
passos, em que n é a ordem do sistema de equações. 

Convém, antes de continuar, determinar a expres- 
sãodo gradiante g”. 

Derivando (1.2) em ordem a q e particularizando 
para q = q, obtém-se o 


(2.9) 


Tal é a expressão do gradiantea, 


[*) À literatura de língua inglesa chama-lhe usteepest déescent method». 
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3. METODO DOS GRADIANTES CONJUGADOS 


O método diz-se dos gradiantes conjugados porque 
os vectores D' se fazem K-conjugados, isto é, tais que 


Tm 


DIK D=0 parart+s (3.1) 
Vai mostrar-se que, se esta condição se verificar, 
a solução exacta se atinge em n passos, 
Observe-se primeiro que, em virtude de (2.1), 


q=.+ 53 kxD (3.2) 


em que r pode tomar qualquer valor inteiro entre zero 
an—71, 
Multiplicando à esquerda por K, obtém-se 


Kq=koaq+- * KR D (3.3) 


e, oprtanto, em virtude de (2.9), 
A q 
g=g+ 2Z à» kKD (3.4) 


Tira-se de (3.4), considerando a simetria de K, 


n—3 
Gir =94D'+ Z ha Da! k D' (3.5) 


==T 


Impondo agora (3.1), vem 


g!D' =g/D +NxNDTK D' 13.6) 


mt 


e, portanto, em virtude de (2.9), 
ag! D' = O (3.7) 


em que r, lembra-se, pode tomar qualquer valor desde 
zero até n — 1 
Admitindo que os sucessivos vectores D' são linear- 


mente independentes, (3.7) implica 9", vector de n 
alementos, seja ortogonal a n vectores linearmente inde- 
pendentes e, portanto, que 


g'=Q (3.8) 

Em virude de (2.9), tem-se então 
K q =0 (3.9) 
O que demonstra o que se pretendia, isto é, que q" 
raiz ce (1.1) e, portanto, a convergência em n passos. 
Resta um problema: como formar uma sucessão 


de vectores D' que sejam K-conjugados? 
Em geral, toma-se . 
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D'= —- gº (3.10) 


como se se tratasse do método da descida mais ingreme. 
Adopta-se seguidamente a fórmula de recorrência 


prri==g7+"+go D' 2.14) 


a determinam-se os 5' de modo que (3.1) se verifique. 
Observa-se que basta impor que cada vector seja 

K-conjugado do anterior para que o seja, no caso de se 

verificar a independência linear, de todos os anteriores, 
impõe-se portanto 


DTTK D=-=-gTIKD+gS DD kODO (3.12) 


o que fornece 


g+HTIK D' 
B' = E TIA (3.13) 
4. RESUMO DE FÓRMULAS 
Parte-se de uma solução inicial q”. 
Tem-se, em virtude de (2,9), 
qe =K - a (4.1) 


im mm 


Determina-se seguidamente Dº, por (3.10), e Nº por 
(231. E 

Obtém-se q' por intermédio de (2.1). 

De novo se determina o gradiante, g', usando 
(2.9), depois A usando (3.11), x! usando (2.7), q 


usando (2.1), e assim sucessivamente, E 

As operações estão esquematizadas no Organo- 
grama |. 

Observa-se que, no caso de equações de sistema 
não serem lineares, a solução exacta não se atinge 
exactamente em n passos, embora, após os primeiros 
passos, a convergência se dê muito rapidamente. O Or- 
ganograma | refere-se a essa situação, prevendo-se por- 
tanto a paragem do processo somente quando o quadra- 
do módulo do gradiante, ou seja, a soma dos quadrados 
dos residuos, seja inferiores a um número : que se 
supõe dado. 


5. APLICAÇÃO A ANÁLISE DE ESTRUTURAS 


A aplicação do método à análise de estruturas não 
ofereceria dificuldades especiais se o utilizássemos de 
modo convencional, isto é, depois de formada a matriz K, 

Mas é justamente a formação da matriz K que se 
deseja evitar. Mais precisamente, deseja-se resolver O 
sistema a partir das matrizes de rigidez elementares, kt, 

Na discussão que vai seguir-se, supõe-se que as 
matrizes kº estão formadas, que os nós da estrutura fo- 
ram dividamente numerados e que faz parte dos dados 


o conjunto das listas dos nós dos vários elementos, isto 
é, uma matriz m de números inteiros, com tantas linhas 
quantos os elementos e com número de colunas igual ao 
número de nós de um elemento. O elemento genérico 
de m, m,;, Feprosesta o número de ordem do nó j do 
elemento e. 

Representando por q* o vector dos deslocamentos 


associado com o nó j do elemento e, por q, o vector dos 
deslocamentos associado com o nó | da estrutura, pode 
eZcrever-se 
q; Um (5.1) 
O vector das forças generalizadas relativas ao nó 
| do elemento e, Q*, pode pois calcular-se por 


(5.2) 


em que Ke j designa uma submatriz genérica da matriz ke. 
O vector das forças generalizadas globais que re- 
sultam dos vectores das forças elementares QO* pode ser 


determinado espalhando e acumulando os vectores Q* 
sobre o vector OQ. Trata-se de uma operação facilmente 


programável, semelhante à que é habitual realizar-se 
para kConstruir a matriz K a partir das matrizes kº. 


im am 


E 

Representa-se essa operação por R, em que E designa 
=| 

o número de elementos, escrevendo-se 


: 
O= RLlQº) (5.3) 


0 


Se em vez do vector Q, dado por (1.1), se perten- 
der obter o vector a 


H=KD' (5.4) 


cm 


= 


poderá proceder-se do mesmo modo, isto é, calculando 
primeiro os vectores 


Hr, = > ke «Dr 
y 


(5.5) 


e em seguida o vector 


H = R (Hr) 


= 


(5.6) 


A análise da estrutura faz-se pois de acordo com 
o Organograma ll, 

Faz-se notar que as matrizes kº pocerão ser memo- 
rizadas em fita magnética ou discos, sendo chamadas à 
memória central uma vez por cada iteração. Nos casos 
de não linearidade fisica as matrizes k” deverão ser 
calculadas de novo em cada iteração, o que não com- 
plicará excessivamente o programa. 


6. CONCLUSÕES 


Corclui-se que é possivel analisar estruturas sem 
chegar a formar a matriz de rigidez global da estrutura, 
isto é, directamente a partir das matrizes de rigidez dos 
elementos. Observa-se que em última análise, o espa- 
lhamento e acumulação sobre uma matriz são substituíi- 
dos pelo espalhamento e acumulação sobre um vector. 

Trata-se de uma conclusão que modifica completa- 
mente os termos em que o problema da programação da 
análise de estruturas costuma ser posto. 

À experiência dirá se o novo método oferece reais 
vantagens, tanto nos casos lineares como nos não-linea- 
res para os quais parece muito especialmente adequado. 
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UM MÉTODO DE FORÇAS FICTÍCIAS 
PARA A RESOLUÇÃO DE PROBLEMAS 
GEOMÉTRICAMENTE NÃO LINEARES 


RESUMO 


Apresenta-se um método para a resolução de pro- 
blemas geometricamente não lineares caracterizado pela 
utilização de forças fictícias obtidas sem qualquer recur- 
so a integrações sobre os elementos ou à forma não- 
“linear das equações deformações-deslocamentos. 


1. INTRODUÇÃO 


O objectivo da presente comunicação é a apre- 
sentação de um método para a resolução de problemas 
geometricamente não-lineares, o qual parece ter vanta- 
gens substanciais em relação a outros de que os auto- 
res têm conhecimento. 

Na resolução de problemas geometricamente não- 
“lineares há dois aspectos a considerar. O primeiro é o 
da formulação do problema, isto é, o da escrita das 
equações que têm de ser resolvidas. O segundo é o 
da resolução dessas equações, 

À presente comunicação refere-se essencialmente 
ao primeiro aspecto. Quanto ao segundo, que só nos 
interessou passageiramente ao elaborarmos os exemplos 
numéricos, deixamo-lo para posterior consideração, 
Observamos contudo que a experiência colhida com ou- 
tros métodos pode ser facilmente aproveitada e que o 
assunto tem sido objecto das atenções de vários auto- 
res, entre os quais mencionamos Oden [1] e Mar- 
çal. [1] 

A apreciação das várias formulações que têm ser- 
vido de base à resolução do problema numérico está 
também fora do âmbito do presente trabalho. 

Indica-se no entanto que uma revisão e sintese dos 
vários métodos que têm sido utilizados foi recente- 
menta apresentada por Stricklin e outros. [2] 

Na comunicação de Stricklin são essencialmente 
considerados dois tipos de formulações: a formulação 
totalmente Lagrangeana, que utiliza as coordenadas dos 
pontos na configuração inicial, e o método das coorde- 
nadas móveis, em que o sistema de coordenadas vai 
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SUMARY 


A method for the solution of geometrically non-li- 
near problems is presented which is characterized by the 
use of fictitious forces obtained without resorting to inte- 
grations over the elements or to the non-linear form of 
the strain-displacement equations. 


variando com a carga, e que só não é denominado 
formulação Euleriana porque a variação da carga se faz 
de modo discreto, por incrementos finitos, dentro de 
cada incremento, utilizada a formulação Lagrangeana. 

No primeiro caso, a consideração dos termos não- 
-lineares nas relações deformações-deslocamentos que, 
no caso de Elasticidade tridimensional, são 


> (uu Ui) (1) 
conduz à utilização de forças fictícias, funções dos des- 
locamentos, cujas expressões envolvem integrais esten- 
didos ao domínio. 

No segundo caso, a consideração dos mesmos ter- 
mos conduz a uma matriz de rigidez dita geométrica, 
que se adiciona à matriz de rigidez elástica e cujos ele- 
mentos são também calculados por meio de inte- 
grações. 

O método que seguidamente se apresenta pode 
considerar-se Euleriano. Conduz também à introdução 
de forças fictícias, mas o cálculo destas é muito mais 
simples que na formulação totalmente Lagrangeana já 
referida. 


2. DETERMINAÇÃO DAS FORÇAS FICTÍCIAS. CASO 
GERAL 


Seguindo a notação utilizada noutros trabalhos do 
autor, [3] as equações de qualquer modelo estrutural 
podem tomar a forma 


equações de equilibrio ............... o E s=T 42) 


equações deformações-deslocamento... 


cd 

|| 
ur] 

E 
to 
gi 


f 


equações tensões-deformações ......... s=H a (4) 


desde que se admita, para simplificar, a linearidade 
fisica. 

E e D são operadores diferenciais e H é uma ma- 
triz simétrica de coeficientes. Denotam-se por s, 8 6 u os 
vectores das tensões, deformações e deslocamentos. 

Se as equações forem escritas, como se supõe 
(formulação Eureliana), para um sistemade coordenadas 
associado ao corpo deformado, o operador E é linear. 


Quanto ao operador D, só poderá transformar-se num 
operador linear D, se se desprezarem certos termos, 
Isto pode fazer-se, no caso de tanto as deformações 
como as rotações serem muito pequenas, isto é, no caso 


de ser válida a chamada hipótese dos pequenos deslo- 
camentos. 


É sempre possível, no entanto,, decompor o campo 
u em duas parcelas, U e à, 


1] = 


| E 


+ à (5) 


das quais WU corresponde a pequenas deformações e 
rotações, Nessas condições, tem-se 


Du=D 0+D (6) 


- ad ia 


Associando (3), (5) e (6), obtém-se pois 


e=D, (lu—W) +D à (7) 

Eliminando as tensões e deformações no sistema 

da equações (2), (4) e (7), conseguem-se as equações 
de Navier modificadas 


E HD u=t+f, (8) 


ima 


Bm que 
H=-EM D-DuA o 


é o vector das forças fictícias. 

Claro que nos dificulta a utilização oa equação (9) 
é a determinação do campo à, isto é, de um campo 
de deslocamentos suficientemente próximo de u para 
que a diferença seja um campo de pequenos desloca- 
mentos. 

O problema resolve-se se se considerar um suces- 
são de estados, suficientemente próximos uns dos outros 
para que, em cada um, possa tomar-se para campo àÔ 
o campo de deslocamentos do estado anterior, sendo 
conhecido o do estado inicial. 

E o que acontece se se supuser a carga variando 
por sucessivos incrementos a que correspondam peque- 
nas variações dos deslocamentos. O campo de desloca- 
mento à inicial pode então tornar-se como nulo, se nula 
for também a carga inicial. 


Observa-se que, procedendo deste modo, o pro- 
blema não-linear se transforma numa sucessão de pro- 
blemas lineares, visto que o operador D,, que apa- 
rece no primeiro membro da equação (8) é um opera- 
dor linear. 

O problema simplifica-se ainda no caso especial 
de o campo de deslocamentos à ser um deslocamento 
dae corpo rígido, o que poderá acontecer se as deforma- 
ções forem muito pequenas e as rotações tiverem pe- 
quenas variações no domínio A correspondente ao corpo, 
isto é, se o campo de deslocamentos do corpo estiver 
muito próximo de um deslocamento de corpo rigido. 

Então, com efeito, 


D à=0 (10) 


Ei 


e a expansão das forças fictícias é 
F=-E H Dô (11) 


A simplificação provém sobretudo de as forças fic- 
tícias poderem ser calculadas sem que seja necessário 
conhecer ou utilizar as relações deformações-desloca- 
mentos não-lineares. Basta conhecer o operador D, , isto 
é, a forma linear das relações deformações-desloca- 
mentos, 

O campo ú poderá ainda ser definido seccional- 
mente, subdividindo o domínio em sub-domínios sufi- 
cientemente pequenos para que, dentro de cada um 
deles, o campo de deslocamentos esteja muito próximo 
de um deslocamento de corpo rigido, desde que se 
tome evidentemente o cuidado de não prejudicar a con- 
tinuíidade das componentes do deslocamento total sobre 
as fronteiras dos subdominios. 

Este último processo torna-se especialmente van- 
tajoso se se discretizar a estrutura, ernpregando o mé- 
todo dos elementos finitos. 


3. APLICAÇÃO AO MODELO DISCRETO 


No modelo discreto, o domínio é um conjunto de 
subdomínios correspondentes a elementos finitos. 

Para cada elemento, podem considerar-se vectores 
das tensões, deformações e deslocamentos, s”, e“' e uº 


e um vector das forças exteriores actuantes no elemento, 
fº. As equações válidas para um elemento genérico, 


escritas para o estado de deformações que se pretende 
analisar, tomam a forma 


equações de equilibrio ..... Di pa 1 E sº=f* (12) 
equações deformações-deslocamentos e'=D' u” (13) 


equações tensões-deformações ...... s"=Hº e" (14) 


Admitindo a hipótese dos pequenos deslocamentos, 
pode tomar-se, em vez de D*, um operador linear D' | 
É então possível demonstrar que E 


E = Dº,! (15) 
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Às equações (12), (13) e (14) há que associar as 
equações de compatibilidade 


=" (16) 
a de equilibrio 
Mo pPmef (17) 
E DE 
globais da estrutura. 


Não sendo válida a hipótese dos pequenos deslo- 
camentos, u' pode decompor-se em iduas parcelas, 


ue 0º, 
W=0" + (18) 
representando úº um campo de pequenos deslocamentos. 
Tem-se | 
D'u =D à + Dºd (19) 
e, portanto, 
e Doo qu — Db) =+ Dº Õº (20) 


Combinando (12), (14) e (20), e atendendo a (15), 
obtém-se 


E ="+E rp nº Dra (21) 


mão 


em que 


k* = Dr! Nº Dº (22) 
é a matriz rigidez do elemento. 

A equação (21) pode simplificar-se se as defor- 
mações forem muito pequenas e se as dimensões do 
elemento forem suficientemente limitadas para que, 
dentro dele, as variações das rotações sejam também 
muito pequenas. 

Então, com efeito, o deslocamento uº estará muito 
próximo de um deslocamento rigido que poderá ser 
tomado para deslocamento 0º. Tem-se pois 


asa 


D*'0=0 (23) 


a 


a a equação (21) simplifica-se para 


erPr=P+HRi O (24) 


ca a ni 


introduzindo (16) em (24) e (24) em (17), obtém- 
-58 finalmente 


ZA KA) U=F+HF, (25) 
O e e a é 


sendo 


FDA td (26) 


et Dn im 


o vector das forças fictícias. 
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Este vector pode pois ser determinado calculando 
primeiro, para cada elemento, o vector 


Pe=k O (27) 
das forças fictícias elementares associadas com o des- 
locamento rigido do elemento, à", e, seguidamente, de- 


terminando para cada nó da estrutura as resultantes 
das forças nodais elementares dos nós dos elementos 
que com ele coincidam, exactamente como se se tra- 
tasse de forças reais. 


No que se refre ao cálculo do vector 6º, há que 


repetir o que se disse na Secção anterior. Poderão pois 
considerar-se incrementos de carga, definindo vários 
estados sucessivos e, em cada um destes, tomar para 


deslocamento de corpo rigido, 0º, o definido pela trans- 


lação e rotação aproximadamente médias do elemento 
no estado anterior. 


4. PROGRAMA DE CÁLCULO AUTOMÁTICO 


Com o objectivo de experimentar o método ex- 
posto em 3, elaborou-se um programa de cálculo auto- 
mático que se designou por METFORFIC (método das 
forças fictícias). 

Numa primeira versão, O programa permite a aná- 
lisas geometricamente não linear (grandes deslocamentos 
mas pequenas dsformações) de estruturas articuladas 
planas com barras rectas de secção constante, Numa 
segunda versão, as barrras passaram a apresentar rigidez 
de flexão. 


Supõe-se que a carga varia por incrementos man- 
tendo-se, em cada estado, proporcional (factor de pro- 
porcionalidade «) a um vector de forças exteriores F., 


ia 


Dos dados fazem parte: 


— número total de nós, apoios e barras; 

— módulo de elasticidade; 

— coordenadas de nós; 

— número dus nós extremos de cada uma das 
barras; 

— comprimentos das barras e áreas das secções 
transversais respectivas; 

— componentes do vector F; 


— lista dos valores de a; a seguir a cada um é 
necessário fornecer um dos valores boleanos 
T ou F, conforme se deseja que sejam ou não 
impressos os resultados correspondentes. 


A lista de dados apresentados presupõe evidente- 
mente uma numeração prévia dos nós da estrutura. 
impõe-se a condição de que os apoios sejam numerados 
em último lugar. 

Observa-se que os valores de w podem ser repe- 
tidos. A tais repetições correspondem iterações que 
podem melhorar a solução correspondente ao valor de 
a respectivo, sendo possivel impor, ou o número de re- 
petições ou o erro máximo admissível. 


5. EXEMPLOS DE APLICAÇÃO 


5.1 Determinar o diagrama F - à relativo à estrutura 
representadas na fig. 1 


A estrutura é constituída por duas barras articula- 
das nas extremidades. À força F é vertical e actua no 
ponto de cruzamento das duas barras. O deslocamento 
5 é o deslocamento do ponto de aplicação da força. 


[)= 3 | E 
DR 

+] = Di VR a qu tm o hm aii e 
e ee = | 10º ts pm - 
3] | E 

| RE 

100 em 
Fig. 1 


Trata-se de um problema cuja solução exacta se 
determina facilmente, pelo que é tradicional escolhê-la 
para o ensaio de técnicas de tratamento de problemas 
geometricamente não-lineares. 

Tendo em vista que a um mesmo valor de F podem 
corresponder vários valores de 5, e que o diagrama F-à 
apresenta pontos em que a tangente é paralela ao eixo 
dos é, não convém fixar incrementos da força mas 
incrementos do deslocamento, 

Como o programa METFORFIC não permite a fi- 
xação de incrementos de deslocamentos, recorreu-se ao 


—— SOLUÇÃO EXACTA 
- METFORFIC 


-+"000="— + FORÇAS FICTÍCIAS NULAS 


artifício de considerar uma barra suplementar vertical, 
da rigidez muito grande em comparação com a da estru- 
tura inicial (fig. 2). 


NGS h] 
UA 


t ro 


Fig. 2 


A força total F, passou a ser absorvida pela barra 
suplementar (parcela F.) e pela estrutura inicial (par- 
cela F). Como a barra é de muito maior rigidez, F é 
praticamente igual a Fr, eo deslocamento 5 é prati- 

Fr ', 
camente igual a (E) | em que (E9)) e |, são respecti- 
"8 
vamente, a rigidez e o comprimento da barra suple- 
mentar. Impor F, é pois praticamente impor o desloca- 
mento, e o diagrama F-ô da estrutura inicial (fig. 1) 
é o diagrama F-ã da estrutura transformada (fig. 2). 

Procurou-se conseguir em &B0 passos o desloca- 
mento máximo de 400 cm, obtendo-se o diagrama da 
fig. 3. 

Na fig. 4 representa-se a variação da força fictícia 
com o deslocamento. Pode observar-se que a força fictícia 
toma valores muito maiores que a força F efectivamente 
aplicada. 


Fig. 3 


10 
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-— I)UD0-— 


59000 


plsjagidao, 


Fig. 4 


Apresentam-se ainda na fig. 3 os resultados (cru- 
zes) que se obteriam se, embora considerando a pro- 
gressiva alteração de geometria da estrutura, para efei- 
tos de cálculo das sucessivas matrizes de rigidez, se 
supusessem nulas as forças fictícias. 


5.2 Determinar a deformada de uma barra de sec- 
ção constante encastrada numa extremidade e 
solicitada na outra por um momento (fig. 5). 


Trata-se também de um problema de solução bem 
conhecida. Sabe-se que a teoria linear atribui à consola 
a forma de uma parábola, e que a forma exacta, dentro 
do quadro da Teoria da Elasticidade com deslocamentos 
finitos, é a circular. 


MH =º Rg em 


Q=W em? 
Ll= 5 emt 


E = Tx” K9,/4m? 


Fig. 5 


Dividiu-se a consola em 10 elementos e atingiu-se 
o momento aplicado de 10º Kgcm, correspondente à 
rotação de 1 radiano do seu ponto de aplicação, por 
meio de 10 incrementos iguais da carga. Em cada nível 
de carga efectuaram-se iterações até se reduzir o erro 
relativo (*) dos deslocamentos a 10— 4.0 número de 
derá ter interesse muito especial no estudo de estados 
iterações consta do Quadro |: 


QUADRO | 
N.º de ordem do incremento | 1 | 2 | 3 | 4 |8 | | 7 | 8 | 9 | 10 
N.º de iterações | 3 | 4 | 4 | 5 | 5 | 5 | 6 | [) | 7 | 8 
Obtiveram-se para valores finais das coordenadas | ML? ML 
da extremidade da consola os valores que, no Quadro Il, linear (u 2E| à VEM | = E] ) 


são comparados com os exactos e o obtidos pela teoria 


(*) Isto é: fez-se [u 
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7 | & | u” | 10 ”*, sendo u!? o deslocamento correspondente à iteração de ordem m. 
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QUADRO II 


Valores obtidos 


Valores exactos 


Teoria linear 50 


y (em) a 
84,1 8 E 0,99999993 
84,15 1,00 E 
100 ao = 


As trajectórias dos nós dos elementos e as defor- 
madas sucessivas podem ser apreciadas na fig. 6. Com- 
pare-se a deformada final com a correspondente à teo- 
ria linear. 


ex (10º Kg) 


Fig. 7 


6. CONCLUSÕES 


O presente trabalho pretende unicamente apresentar 
o método e, com o auxílio de dois exemplos simples, 
demonstrar a sua aplicabilidade. 

Requere-se agora uma exploração sistemática das 
respectivas possibilidades, particularmente no trata- 
mento de problemas de instabilidade. Julga-se que po- 
derá ter interesse muito especial no estudo de estados 
de pós-encurvadura. Impõe-se também uma comparação 
com outros métodos até agora utilizados. Tudo isso 
fica porém para trabalhos futuros, 
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M (10º Kgem) 


A fig. 7 mostra como as forças fictícias evoluem 


coma carga. 


Fe. (10º Kg) 
8) fy 


mg 


== 


o 
É 
-— 10 E 


Me (10º Kgem) 


M 0º Kg em) 
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INTRODUÇÃO À MECÂNICA DA FRACTURA 


RESUMO 


A mecânica da fractura relaciona as propriedades, 
defeitos de compacidade dos materiais e suas condições 
de aplicação como o fenómeno da fractura. 

Neste trabalho apresentam-se conceitos e formula- 
ções analíticas básicos da mecânica da fractura. 


1. TERMINOLOGIA BÁSICA E OBJECTIVO 


Com o fim de minimizar ambiguidades conside- 
rou-se conveniente precisar o significado dos termos 
utilizados. 

Utilizaremos a palavra fractura a propósito de casos 
em que a ruptura se dá essencialmente devido à propa- 
gação duma fenda, característica que as distingue dos 
putros casos em que a ruptura se dá essencialmente 
devido a excessiva deformação plástica. 


Designaremos por mecânica de fractura a ciência 
que trata da quantificação dos factores que condicionam 
a extensão duma fenda em função das tensões e defor- 
mações. Ao aprofundar-se tal ciência depara-se rapida- 
mente com a necessidade da utilização duma grandeza 
correspondente à resistência que o material opõe à ex- 
tansão duma fenda, grandeza essa que designaremos por 
tenacidade à fractura ou simplesmente tenacidade, 


A mecânica da fractura não se limita ao tratamento 
da extensão instável de fendas provocada simplesmente 
palas forças exteriores, mas pode igualmente aplicar-se 
à extensão de fendas provocada por fadiga, corrosão 
sob a tensão, fluência, fragilização pelo hidrogénio ou 
por combinação de dois ou mais destes mecanismos. 


A aplicação dos processos analíticos da mecânica 
da fractura ao projecto duma obra com o objectivo de 
garantir a segurança dos seus componentes contra a 
ruptura designa-se por controle da fractura. O control 
da fractura implica nomeadamente a especificação da 
tenacidade do material, dos processos de fabricação 
e dos defeitos admissíveis. 

O controle da fractura propõe também indicar para 
cada elemento duma obra da material, dimensões a gao- 
metria definidas, quais os defeitos admissíveis para 
que o dito elemento resista a condições de serviço bem 
determinadas. 
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ABSTRACT 


The fracture mechanic relates the properties, the 
compacity faults of materials and its condition of apli- 
cation with the fracture phenomenon, 

This work presents basic concepts and analytic 
formulations of the fracture mechanics. 


A luz da mecânica da fractura a determinação das 
dimensões dos defeitos é feita por via analítica entrando 
em conta com factores geométricos, com as tensões 
de serviço e seu tipo, e com parâmetros característicos 
do material nas condições de serviço e espessura pre- 
sentes, nomeadamente com a tenacidade, À tenacidade 
passa assim a ser um parâmetro que se pode e deve 
utilizar no cálculo, logo na fase de projecto da obra, 
tal como por exemplo a tensão limite de elastigidade. 
Esta situação contrasta com a tradicional em que a tena- 
cidade, ainda vvulgarmente expressa em resiliência Char- 
py, só pode ser utilizada duma forma empírica, por 
comparação com a de obras do mesmo tipo, material e 
condições de utilização que prestaram boas provas em 
serviço. 

Compreende-se que a situação tradicional possa 
conduzir quer à falta de segurança quer à falta de eco- 
nomia quer á falta de segurança. À falta de economia 
pode surgir quando sendo a comparação atrás referida 
possível, se adoptaram critérios mais drásticos do que 
na realidade seria necessário. À falta de segurança pode 
surgir quando se trata de obras que impliquem uma 
combinação materiais-espessura-condições de serviço 
tais que o projectista não possua dados estatísticos so- 
bre obras similares; é o que pode suceder quando se 
passa para materiais mais resistentes, ou para maiores 
espessuras, ou para utilização a mais baixas tempera- 
turas, etc. Não é assim de admirar o interesse crescente 
que por toda a parte está despertando o estudo da me- 
cânica da fractura, Para este interesse contribuiram tam- 
bém os problemas encontrados nos projectos aeroespa- 
ciais com os materiais de alta resistência e os inerentes 
à utilização de fortes espessuras nos componentes de 
centrais nucleares. 

A mecânica da fractura conta já hoje com êxitos 
na aplicação quer a construção de materiais de elevada 
resistência e/ou forte espessura, quer a construções de 
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materiais e espessuras mais correntes como sejam reser- 
vatórios sob pressão construidos de aços ao C-Mn de 
média resistência. 

Note-se que as dimensões dos defeitos admissíveis 
podem ser, consoante os casos, muito diversas, desde 
dezenas de centimetros até dimensões inferiores ao mi. 
nimo detectável pelas actuais técnicas de controle não 
destrutivo, 

Neste último caso, correspondente a situações ex- 
cepcionais, só o ensaio de sobrecarga poderá dar uma 
ideia da capacidade da construção resistir aos esforços; 
mas então a mecânica da fractura de nada serviu — 
para ser utilizada nestes casos excepcionais é nacessá- 
rio que se desenvolvam as técnicas de control não 
destrutivo. 


2. CARACTERÍSTICAS GERAIS DAS FRACTURAS 


O fenómeno da fractura pode processar-se por me- 
canismos de clivagem, quase-clivagem, coalescência de 
microcavidades, escorregamento a estricção, fadiga, cor- 
rosão sob tensão, fragilização pelo hidrogénio, fragili- 
zação devida ao contacto com outros metais líquidos, 
ou por combinação de dois ou mais destes mecanismos, 
que é de resto O caso mais frequente. 

A escala macroscópica as superfícios resultantes 
da fractura, apresentam essencialmente três tipos dis- 
tintos: 


a) Plano — quando a superficie resultante da frac- 
tura é plana e normal à tensão que a originou; 


b) inclinado — quando a superfície resultante da 
fractura forma um ângulo de cerca de 45º com 
direcção da tensão que a originou; 


c) fibroso — quando a superficie resultante da 
fractura apresenta o aspecto de fibras, prove- 
nientes de estricção acentuada do material, 
orientadas segundo a direcção da tensão que as 
originou. 


Mais frequentemente em vez destes tipos puros 
ancontra-se o tipo mixto, correspondeanta à combinação 
daqueles. 

O tipo da superfícia depende dos mecanismos mi- 
croscópicos de fractura intervenientes, Por sua vez fac- 
tores macroscópicos determinam os mecanismos micros- 
cópicos. Ássim por exemplo o aumento da espessura em 
provetes de aço pode fazer com que a fractura por cli- 
vagem predomina sobre a fractura por coalescência de 
mitrocavidades; o mesmo se passa quando, por exem- 
plo em provetes de titânio, a acuidade do entalhe au- 
menta, 

As fracturas planas podem ser geradas por quais- 
quer dos mecanismos citados. As fracturas inclinadas 
são quase sempre geradas pela coalescência de micro- 
cavidades à qual se segue a ruptura por escorregamento 
e estricção; todavia a ruptura também se pode consumar 
por fadiga ou por clivagem; este último caso sucede 
por exemplo em provetes de Cr e de grão fino que após 
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astricção severa sa rompem por clivagem escorvada 
em diversos pontos da secção reduzida por estricção. 

As estruturas mais clássicas da mecânica da 
fractura não se baseiam nos mecanismos microscópicos 
de fractura citados mas sim num mecanismo macroscó- 
pico. Falar dum mecanismo macroscópico de fractura 
é aparentemente uma abstracção, visto que os defeitos 
macroscópicos são resultantes dos processos microscó- 
picos implicados na fractura. Já a teoria da elasticidade 
assenta numa abstracção ao considerar os corpos elás- 
ticos como homogéneos e isótropos. 

Segundo a mecânica da fractura mais tlássica, 
uma fenda propaga-se quando os seus lábios ao afas- 
tarem-se criam, na região do fundo da fenda, concentra- 
cos da tensões e de extensões que o material é incapaz 
de suportar. Avalia-se assim a situação do material no 
fundo da fenda através das grandezas macroscópicas 
(tensões e extensões), Mas estas grandezas macroscópi- 
cas reflectem a situação microscópica, Compreende-se 
assim Oo êxito conseguido pela aplicação duma teoria 
baseada em hipóteses macroscópicas. 


3. ESTADOS PLANOS DE TENSÕES E ESTADOS PLA- 
NOS DE EXTENSÕES 


Na continuação referir-nos-emos a estados planos 
de tensões bem como a estados planos de extensões. 
É pois conveniente rever as características de cada um 
daqueles estados, 

Reportar-nos-emos desde já ao caso que mais nos 
interessa isto é de chapas e a um sistema de eixos x, 
Y. Z, sendo o eixo dos zz normal à superfícia da chapa. 

Diz-se que se está em presença dum estado plano 
de tensões quando se verifica 5, — Ta — Tay = 0 

Este estado ocorre por exemplo quando se submete 
uma chapa muito fina a tensões uniformes de tracção 
orientadas segundo o plano da chapa. Na verdade qual- 
quer «fatia» paralela às faces da chapa se pode consi- 
derar livre de tensões, nas suas superfícies normais ao 
eixo dos zz por estas se encontrarem muito próximas 
das faces da chapa que por hipóteses estão livres de 
tansõas. 

Diz-sa que se está em presença dum estado plano 
de extensões quando se verifica :, — Tx Tay” O 

Este estado ocorre por exemplo na zona central 
duma chapa muito espessa submetida a tensões unifor- 
mes de tracção orientadas segundo o plano da chapa, 
ao passo que as zonas laterais, nas proximidades das 
superfícies da chapa, ficam sujeitas a um estado plano 
de tensões. Na verdade à variação da espessura da zona 
central opõem-sa as zonas adjacentes; o mesmo não 
suceda com as zonas laterais às quais se aplicam racio- 
cínios idênticos aos feitos sobre as chapas finas. 

Em provetes prismáticos entalhados num ou em dois 
lados, pode encontrar-se, após fractura por aplicação 
de esforços de tracção uniformes e normais ao entalhe 
uma fractura do tipo mixto, plana na zona central do 
provete a inclinada nas zonas laterais. A forta triaxili- 
dade existente na zona central impediu a actuação das 
tansões da corte dando-se a fractura porque a tensão 
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principal de tracção, normal ao plano do entalhe, ul- 
trapassou um valor crítico; sob a acção desta tensão a 
fractura propaga-se segundo o plano do entalhe. Nas zo- 
nas laterais pelo contrário a reduzida triaxilidade per- 
mitiu a actuação das tensões de corte o que explica o 
tipo inclinado de fractura. Estamos pois em presença 
do caso já citado em que na zona central existe 
um estado plano de extensões e nas zonas laterais 
estados planos de tensões. 

Sendo assim e supondo que o plano x z, coincide 
com o plano que contém o entalhe ter-se-á na zona 
central: 


donda 


o =v(8, + 0.) 

No plano de entalhe tem-se pela teoria da elasti- 
cidade (válida para pontos situados fora da pequena 
zona plástica que se forma na extremidade do entalhe) 
CE Oy de forma que as tensões principais são 


É ==) 
Y x 


o. = 241, 


4. ANÁLISE DAS TENSÕES E DAS DEFORMAÇÕES 
NA VIZINHANÇA DA EXTREMIDADE DUMA FENDA 


À teoria da elasticidade permite estudar os campos 
de tensões e de deformações existentes na vizinhança 
da extremidade duma fenda, Assim Irwin demonstrou 
que em certas hipóteses de geometria e de solicitações 
as tensões existentes nas regiões das extremidade são, 


em condições elásticas, dadas pelas expressões 
K, 6 | 8 30 
E = Cos — | 1— Bh. 961 ——- 
Px — (Qar) lt 2 2 2 
K, [5 8 39 — 
Ou — — 808 - 1+5sen — sen 
y (2mr)!'? 2 


K | 8 [E] 3a 


é (2arr)p 1)? 
Ga — * lgy 1 q! em estado plano de extensões 


G — 0 em estado plano de tensões 


sendo r e à as coordenadas polares e K, um parâmetro 
denominado factor de intensidades de tensões, 

Dum modo mais geral pode dizer-se que as ten- 
sões nas regiões das extremidades são dadas por ex- 
pressões do tipo 
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K 


E E (2qr) '? f; (5) 


e os deslocamentos são dados por 


r 


u — K, V (gi) (1) 
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K, depende da geometria, do tipo de solicitações 
a das proprisdades elásticas do material, Todavia em 
certas hipóteses de solicitações K, é independente das 
propriadades elásticas do material. 

Note-se que K, não depende do estado em causa 
ser um estado plano de tensões ou de extensões. Ainda 
não se deverá confundir K, com um factor de concen- 
trações de tensões pois ele depende destas. 


b. EXTENSÃO DUMA FENDA EM SÓLIDOS FRÁGEIS 
5.1 Introdução 


Apresentar-se-ão a seguir diversos critérios de ex- 
tensão de fendas em sólidos frágeis. 

Quando se verificam as condições correspondentes 
ao critério aplicado, a fenda alonga-se isto é, sofre uma 
extensão. 

A alteração da geometria resultante da extensão 
da fenda provoca directamente uma alteração das ten- 
sões reinantes nas regiões das extremidades. Há tam- 
bém uma alteração provocada indirectamente pela ex- 
tensão da fenda, De facto a extensão da fenda permi- 
tindo um aumento das dimensões na direcção normal 
à fenda, pode provocar uma redução das solicitações 
correspondentes: é o que se passa por exemplo num 
ensaio em que é fixa a distância entre as garras da má- 
quina, ou num recipiente contendo um fluido incom- 
pressivel quando, em virtude da extensão da fenda, se 
dá um derramamento de fluido, 

Por outro lado a própria tenacidade pode manter-se 
constante, mas também pode aumentar ou diminuir 
quando o comprimento da fenda aumenta, Uma das 
causas desta alteração reside no facto dos campos de 
tensões e extensões existentes nas regiões das extre- 
midades da fenda variarem durante a extensão desta. 

O que se disse já permite concluir que do facto 
de se verificarem as condições necessárias à extensão 
duma fenda não se pode concluir que se atinja a frac- 
tura. Os critérios de extensão apenas permitem concluir 
que esta se inicia, por isso são vulgarmente chamados 
critérios de iniciação, 

O estudo das fases posteriores à extensão isto é 
a propagação, não será abordada neste trabalho, nem 
para sólidos frágeis nem para sólidos semi-dúcteis e 
dúcteis. 


D.2 Critério de Griffith 


Griffith demonstrou que uma fenda de secção recta 
elíptica achatada de comprimento 2 a, existente numa 
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placa infinita submetida a uma tensão uniforme oc nor- 
mal ao plano da fenda, se propaga de forma brutal quan- 
do o valor da tensão aplicada atinge um valor crítico 
dado por 


em estado pleno de tensões 


=| 
m 


ou por 


2 Er [ii 
gu = | — ROc:*:6 2 
c ct egii E ) 


em estado plano de extensões 


Na dedução destas expressões admitiu-se não ha- 
ver deformação plástica nas extremidades da fenda, 
isto é considerarem-se corpos verdadeiramente frágeis. 
Além desta hipótese não realista a teoria de Griffith 
está sujeita a várias críticas que adiante se apontarão. 
Apesar disso a sua importância para o estudo dos sóli- 
dos frágeis é indiscutível. Por outro lado deu origem à 
mecânica da fractura 


5.3 Critério de Griffith — Orowan 


Reexaminando o critério ce tensão de Griffith, 
Orowan demonstrou com hipóteses idênticas, mas su- 
pondo que o raio de curvatura nas extremidades da 
fenda é da mesma ordem que o raio atóúmico que a 
tensão crítica será dada por 


Eq 


Mais importante porém é a contribuição dada por 
Orowan para a adaptação da teoria de Griffith aos me- 
tais. Nestes a deformação plástica nas extremidades da 
fenda é relativamente importante ao contrário do que 
suceda com os sólidos ditos frágeis. Segundo Orowan 
as expressões de Griffith serão mais correctas se em vez 
do parâmetro y se considerar não só a energia superficial 
T s - mas também a energia necessária à deformação 
plástica Yo (correspondente a um aumento unitário 
do comprimento da fenda) nas extremidades da fenda, 
Demonstra-se que quando o comprimento da fenda é 
macroscópico y, é muito elevado em relação y,. Sendo 
assim a tenacidade dum metal é sobretudo condicionada 
pela sua ductilidade. 

Nos metais a em condições de comportamento 
frágil, y, & Yp estão relacionados com a força de 
extensão da fenda G,e por 


em estado pleno de tensões 


Gic = Ta =| Tp (4), 
sendo Gjç definida como uma força fictícia que pro- 
duza o trabalho correspondente às energias superficial 


a plástica necessárias à expansão unitária da fenda. 
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Ora nas expressões de Griffith e na similar de 
Orowan, aplicáveis aos sólidos ditos frágeis, o trabalho 
correspondente à expansão unitária da fenda intarvém, 
através de y, formalmente reduzido a metade. Então 
para considerar a contribuição do Orowan dever-se-á 
substituir nas expressões citadas y por G,ç /2. 

Então a primeira e a segunda expressão de Griffith 
poder-se-ão excrever respectivamente. 


EGjc 1 
(5) em estado plano de tensões 


: EG ,c 12 
é P (1—»º) a | 


(6) em estado plano de extensões 


Para tensões o inferiores à crítica define-se da for- 
ma análoga G por 


EG, Fê 


- [- | (7) em estado plano de tensões 


“1 


q 
o] 


ES Tm 
| T (1 = 5) a 


(8) em estado plano de extensões 


Atendendo a (4) também se pode definir G je 
como sendo a energia, reduzida à unidade de espessura 
da placa, necessária ao aumento unitário do comprimento 
da fenda. No início da extensão da fenda tal energia 
é investida em trabalho da deformação plástica e em 
energia superficial. 


5.4 Critério da tenasidade crítica 


As tensões principais reinantes na região que en- 
volve as extremidade duma fenda são função da inten- 
sidade de concentração de tensões. Assim para uma 
fenda elíptica de semi-comprimento a e raio de curva- 
tura nas extremidades p, existente numa placa infinita 
submetida a tensões q normais à direcção da fenda, 


demonstra-se que a tensão máxima cms, é dada por 


ú + 2( ) franco + z ( a 


factor de concentração de tensões. À expressão anterior 
lava à conclusão que para fendas muito estreitas, com 
p=>0, O max Seria infinito. Ora na realidade mesmo 
nos metais ditos frágeis dá-se antes da ruptura uma 
certa deformação plástca qua torna finito o raio de cur- 
vatura das extremidades da fenda. 


“máx” 9 um 
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A tensão o mix Não poderá por consequência ser- 
vir de critério de ruptura, 

Talvez pudesse servir de critério a tensão máxima 
atinginda a uma certa distância da extremidade da fenda. 
Imediatamente se levanta o problema de definir esta 
distância. 

É então mais lógico recorrer a um factor definidor 
do campo de tensões na extremidade da fenda. Recor- 
re-se ao factor de intensidade de tensões K, que como 
já se disse depende da geometria, da distribuição das 
forças aplicadas e do comportamento elástico do ma- 
terial. 

O valor de K, aumenta à medida que aumentam 
as tensões aplicadas, Ora quando estas tensões atingem 
um certo valor, que corresponde um valor de factor de 
intensidade de tensões que designaremos por Kje, ini- 
cia-se a extensão da fenda, 

Pode então dizer-se que a extensão da fenda se 
dará quando o factor de intensidade de tensões ultra- 
passar um valor crítico K,c . dependente da geometria, 
da distribuição das forças aplicadas, do material e das 
condições de ensaio. K,c , é pois uma medida da tena- 
cidade do material nas condições de ensaio conside- 
radas. 

Note-se que considerando a mesma geometria 6 a 
mesma distriblição das forças aplicadas, por exemplo 
a distribuição uniforme, K,c dependerá apenas do ma- 
terial e das condições de ensaio. 

O valor de K, pode também calcular-se pela passa- 
gem ao limite: 


K =C lim o y%p 
PU" 


sendo C uma constante 


O cálculo de K, por este processo é geralmente 
dificil. 

Para placas infinitas com uma fenda central, K , 
está relacionado com a tensão aplicada uniforme o e 
com o semi-comprimento da fenda por 


Ko (ga (10) 


e o valor crítico por Kg -— cc y qa 


sd 


(10 a) 


Se as placas de largura finita W ter-se-á para K, 


E q” 
a E 
o (4 )| (11) 


4 


Kj— o 


e para o valor crítico 


ma 4 |! 
Kic = ve] W ol )| 


(11 a) 


Atendendo às expressões (7), (10) e às expressões 
(5), (10 a) obtém-se respectivamente 


C, = ER, (12) 


Kjc = EG,e (12 a) 

ambas em estado plano da tensões, Para o estado plano 
de extensões obter-se-á de modo análogo à custa das 
expressões (8), (10) e (6), (10) 


E 
==) (13) 


1 — 


E 
o 


Ki ri pica a 
IC E = 


(13 a) 


E 


Tal como foi definido, G, é também uma me- 
dida da tenacidade do material, dependendo da geome- 
tria, do material e das condições de ensaio. 


6. EXTENSÃO DUMA FENDA EM SÓLIDOS SEMI- 
-DUCTEIS E DÚCTEIS (*) 


6.1 Análise e Irwin 


Irwin utiliza um modelo que apesar de bastante 
aproximado no que se refere às dimensões da zona 
plástica, é muito importante na prática para a deter- 
minação de K,.. 

As tensões normais a uma fenda achatada de com- 
primento 2 a, existente numa placa infinita submetida a 
tensões uniformes o de tracção, normais ao plano da 
fenda são segundo Irwin, quando se considera que o 
campo de tensões existentes na placa é puramente 
elástico e em estado plano de extensões: 


Gx 
GEE REA (14) 
V C— a” 


Esta expressão pode escrever-se com suficiente 
aproximação sob a forma 


r<<a (15) 


sendo r a distância à extremidade e G um parâmetro re- 
lacionado com a intensidade de tensões nas regiões 
vizinhas das extremidades da fenda, G tem significado 
análogo a G, mas por se referir a sólidos semi-dúcteis 
ou dúcteis que em geral não conduzem a fractura plana, 
omitiu-se o Índice |, 

Os deslocamentos das superfícies das fendas são 
dadas por 


(*) Neste trabalho designam-se por sólidos semi-ducteis aqueles em que no início da extensão da fenda, as zonas vizinhas 
das extremidades daquela se encontram em deformação plástica, embora zonas mais afastadas se encontrem em deformação elástica. 
Designam-se por sólidos dúcteis aqueles em que no início da extensão da fenda se encontram na totalidade em deformação 


plástica. 
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| (às x<a 16) 
7 E V a E us | 


ou de forma aproximada por 


2a 


% = dar = 2 f od r<<a (17) 
E y y q E 
sendo agora r medido no sentico oposto, isto é para 
o interior da fenda. 
As tensões o, estão representadas na fig. 1. Mas 
as tensões para x = a são infinitas o que é impossível. 
Então Irwin, considerando um material elasto-plástico 
puro, admite que dentro da zona de largura R (fig. 2) a 
tensão não ultrapassa a tensão limite de elasticidade o,. 
Admite ainda que fora daquela zona as tensões são as 
correspondentes à distribuição elástica relativa à fenda 
aumentada de r.. Atendendo a que a distribuição de 
tensões no plano da fenda deverá equilibrar a tensão 
aplicada, demonstra-se que R = 2 r;. À distribuição de 
tensões resultante (fig. 2) é pois a correspondente a 
uma fenda fictícia de semi-comprimento a + r,. trun- 
cada no valor Sy Ze: À dimensão da zona plástica 
será segundo (15): 


pi (18) 


6.2 Desenvolvimento de Wells 


Desenvolvendo a análise de Irwin, Wells deduziu 
expressões para o deslocamento da abertura da fenda 
(crack open displacement — €C. O, D.) 6 nas extremi- 
dades, desta para o caso de sólidos semi-dúcteis. 


MODELO ELÁSTICO PURO 
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A fenda ficticia suporta, à distância , das extremi- 
dades, e devido à deformação plástica, um deslocamento 
(aumento) da abertura cujo valor é 5 = 2ny- 

Por outro lado as expressões (15) e (17) podem 
aplicar-se à fenda fictícia desde que se faça r= ,. 
se substitua a por a=a+r, e se faça 6G=G6, 
sendo G, a «força de extensão da fenda» capaz de 
provocar a sua extensão, mas condições consideradas. 
Então substituindo em (17) o valor de r, tirado de (15) 
obtém-se, utilizando as primeiras e as segundas igual- 
dades daquelas expressões, respectivamente: 


fo a 4 G 
35-24, (19) 


Eo q To 


= 


Se consideramos o que se passa desde o início 
de actuação das forças aplicadas, diremos que se criam 
nas extremidades da fenda zonas plásticas que aumen- 
tam progressivamente assim como o deslocamento da 
abertura da fenda. Mas este último não poda aumentar 
indefinidamenta dado que o material tem uma capaci- 
dade de deformação plástica limitada. Dada a forma 
como se definiu G,, o valor de à dado por (19) corres- 
ponde precisamente ao limite isto é ao valor crítico 
a partir do qual a fenda se expande. 

Dentro da zona plástica a energia necessária à 
extensão da fenda engloba a que é absorvida pela de- 
formação plástica e a necessária à criação da super- 
fícia da fenda (energia superficial). Então por um ba: 
lanço energético qualitativo pode escrever-se 


e (19 a) 


A expressão (19 a) é mais geral e por isso prefe- 
rivel a (19) visto que esta última se refere às condi- 
ções elasto-plásticos particulares. 

A partir de (18) e (19 a) pode escrever-se 


| 5 E 5 
Fs = e - 
y Ana, 2 nt id 
a 
O valor de — pode calcular-se a partir de (20) 


F 

y 
ea partir de (15), considerando a forma que esta última 
expressão assume quando aplicada à fenda fictícia 
obtendo-se respectivamente 


(21) 


De (21) conclui-se que os semi-comprimentos da 
fenda e da zona plástica serão iguais quando o = c, 
isto é no início da deformação plástica generalizada. 
Todavia é natural pôr em dúvida que (21) se aplique 
a deformação plástica generalizada. 
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Wells considera que à temperatura de transicção, do 
estado plano de tensões para o estado plano de exten- 
sões, o valor crítico do C. O, D. é dado por 


do = 92m: t (22) 


sendo t a espessura do material, Por outro lado demons- 


trou experimentalmente que quando a temperatura é 
tal que a iniciação da ruptura frágil (*) se dá após, ou 


i 

t 

1 
i 
1 


ZONA PLASTICA 


DISTRIBUIÇÃO APOS 
DEFORMAÇÃO PLÁSTICA 


Fig. 2 


imediatamente antes da cedência generalizada, o valor 
de à correspondente é constante para uma gama extensa 
de valores do comprimento da fenda, quer o ensaio seja 
feito sobre placas largas com uma fenda central quer 
sobre provetes de dobragem com uma fenda lateral. 
Esta constatação não discorda da expressão (22) e de- 
monstra que deve haver alguma relação entre do e as 
características do material provavelmente, como sugere 
Wells, entre do Gata É. 

Wells verificou também por via experimental que 
para temperaturas inferiores à de transicção, corres- 
pondenta a uma mudança do estado plano de tensões 
para o estado plano de extensões: 


a) As dimensões das zonas plásticas existentes nas 
extremiades da fenda tornam-se inferiores à 
espessura da placa. 


b) Os valores do COD correspondentes à fractura 


reduzem-se à medida que a temperatura baixa. 


NOTA 1 


Tal como a G,.. está associado K ., a G,. está 
associado a um factor de intensidade de tensões que 
designaremos por Ka 


Ao passo que K |. se aplica a estados planos de 
extensões, K.. aplica-se a estados planos de tensões. 

O factor de intensidade de tensões dum dado 
material varia em função da espessura tal como se in- 
dica na fig. 3. 


= 
ft (espessura) 


Fig. 3 


Verifica-se que para valores elevados da espessura, 
aos quais correspondem estados planos de extensões 
a fracturas planas, o factor de intensidade de tensões 
se torna constante, sendo assim uma caracteristica do 
material. Nestas condições é designado por Ki 

À medida que a espessura diminui o factor de 
intensidade de tensões aumenta atingindo os maiores 
valores para as menores espessuras correspondentes 
a estados planos de tensões, Nestas condições corres- 
ponde a Ko: O valor K é uma característica dum pro- 
duto. 

O tipo de fractura respectiva é inclinado; para es- 
pessura intermédias, o tipo de superficie resultante da 
fractura é mixto, 


NOTA 2 


O factor K e aplica-se a casos em que a extensão 
da fenda se dá em condições ditas frágeis para as 
quais, as dimensões das zonas plásticas existentes nas 
extremidades da fenda, são diminutas, Ora a dimensão 
da zona plástica é, como as expressões (12 a) e (18) 


o: 
deixam antever, proporcional a | — )º. Por consequên- 


õ 
e 


será aplicável quando a espessura t do ma- 


cia K ãê 


terial for suficientemente superior a = )*: verifica-se 
“a 
Ko 
ser suficiente que t > 2,5 (— )*, sendo então a frac- 
“e 
tura frágil. Se t for inferior mas pouco ao valor fornecido 


pala desigualdade anterior a fractura será semi-dúctil; 
se for bastante inferior será dúctil, 


[*) A expressão «ruptura frágil» está aplicada aqui no sentido macroscópico, significando que a deformação plástica 
medida sobra comprimentos de referência da ordem de grandeza dos habituais é diminuta. 
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6.3 Teoria de Barenbilatt 


A teoria de Griffith apresenta limitações algumas 
das quais se apresentam a seguir, 


a) às deformações na extremidade duma fenda 
saem fora do domínio da elasticidade. Em vir- 
tude da deformação, a elipse por hipótese de 
raio quase nulo nas extremidades, passa a ter 
um raio finito. 


b) O ângulo interno formado pelas faces da fenda 
nas extremidades passa de Oº no início para 
180º depois da aplicação das forças exteriores 
Assim nas regiões que englobam as extremida- 
des da fenda, as partículas vizinhas afastam-se 
muito e dever-se-á considerar o comportamento 
atómico em vez do comportamento dos meios 
contínuos, 


c) No cálculo da redução da energia elástica que 
se verifica quando da extensão da fenda, admi- 
te-se que não existem tensões aplicadas na su- 
perfície da fenda mas por outro lado, admite-se 
que esta superfície possui energia superficial. 
Ora a existência de energia superficial e por- 
tanto de tensão superficial numa superficie 
curva, implica a existência de tensões normais 
à superfície que equilibrem a tensão superfi- 
cial. 


d) A teoria de Griffith considera uma elipse que no 
limita dá uma fenda. Não demonstra que se 
obteriam resultados análogos a partir da passa- 
gem ao limite de outro tipo de descontinuidade. 


Na teoria cujos postulados e conclusões gerais se 
exporá imediatamente a seguir, considerar-se-á indistin- 
tamente a forma inicial que lhe foi dada por Barenblatt 
assim como aperfeiçoamentos posteriores, 

Os resultados obtêm-se pela sobreposição das solu- 
ções dos três problemas de elasticidade seguintes: 


a) Placa sem fenda, e simplesmente sujeita às 
forças exteriores, 


b) A mesma placa, sem forças exteriores, mas 
com fenda, cujas superfícies se encontram su- 
jeitas a uma distribuição de forças iguais 
opostas às que nos locais correspondentes 
criam as forças exteriores na hipótese a). 


c) A mesma placa, sem forças exteriores, com 
fenda cujas superfícies se encontram sujeitas 
às forças de coesão. 

Às forças de coesão (ou de atracção ató- 
mica ou molecular), actuando nas extremidades 
da fenda, tendem a comprimir entre si 
perfícies desta. Consideradas só por si (hipó- 
tese c) induzem singularidades na distribuição 
das tensões de compressão que criam nas extre- 
midades da fenda, ao passo que as forças apli- 


as Sl- 
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cadas, actuando só por si, induzem singularida- 
des na distribuição das tensões de tracção que 
criam nas mesmas zonas, O primeiro postulado 
em que se baseia esta teoria diz que pela sobre- 
posição daquelas distribuições de tensões se 
obtém uma outra não singular. Não havendo sin- 
gularidades na distribuição das tensões tam- 
bém não as haverá nos deslocamentos relativos 
das superfícies da fenda que unirão suavemente 
constituindo um perfil em forma de ponta de 
lança como se indica na fig, 4. 


Fig. 4 


Não há vantagem nesta teoria em considerar que 
as forças aplicadas estão uniformemente distribuidas 
visto se considerar que a distribuição das tensões de 
coesão nas extremidades da fenda não é uniforme. 
Considera-se pois que as tensões aplicadas têm uma dis- 
tribuição mais genérica e impõe-se apenas que ela seja 
simétrica relativamente ao plano da fenda e ao eixo 
normal âquela e que passa pela zona média. 

O segundo postulado em que se baseia a presente 
teoria afirma que as forças de coesão se fazem sentir 
em pequenas zonas a partir das extremidades da fenda 
e para o interior desta. 

O terceiro e último postulado diz que «a forma da 
secção normal da superfície da fenda (em expansão) 
na região das extremidades (e consequentemente a 
respectiva distribuição local de forças de coesão na su- 
perficia da fenda) não depende das forças aplicadas 
e é sempre a mesma para um dado material em condi- 
ções definidas (temperatura, composição e pressão da 
atmosfera envolvente, etc)», 

Segundo esta teoria as forças de coesão começam 
por aumentar à medida que aumentam as forças aplica- 
das até que chegada uma dada altura não conseguem 
mais suportar os efeitos daquelas e a fenda expande-se. 
Mas durante a expansão a distribuição das forças de 
coesão nas extremidades da fenda mantém todas as 
suas caracteristicas assim como a forma da secção 
normal da superficie da fenda nas mesmas zonas — em 
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ponta de lança. Nestas circunstâncias diz-se que a fenda 
se encontra em equilíbrio móvel. Pelo menos no início 
da extensão da fenda não se trata pois duma extensão 
brutal, isto é, rápida mas sim que se dá à medida que 
aumentam as forças aplicadas. 

As bases expostas permitem calcular as forças 
aplicadas capazes de provocar a extensão da fenda a 
partir do conhecimento das forças de coesão, traduzidas 
por uma constanta do material chamada módulo da 
coesão. Permitem igualmente calcular os deslocamentos 
relativos das superfícies da fenda, e por consequência 
a forma desta em cada instante, correspondente a uma 
dada distribuição das forças aplicadas (até a fenda en- 
trar em expansão pois a partir dai a forma das suas 
extremidades é, como já se disse, independente das 
forças aplicadas). 


Note-se que o terceiro postulado, com a introdução 
das forças de coesão distribuidas ao longo de pequenos 
comprimentos nas extremidades da fenda, suplementa 
o modelo continuo da teoria da elasticidade. 


Torna assim possível estabelecer algumas relações 
com os pequenos comprimentos intrínsecos aos modelos 
das redes, isto é com os espaços atómicos ou mole- 
culares. 


O limita que o terceiro postulado impõe às forças 
coesivas torna ainda possivel estabelecer uma relação 
com a resistência teórica do material. 


6.4 Algumas reflexões sobre a teoria de Barenblatt 


Analisando a forma assumida pela extremidades 
duma fenda que se expande de acordo com a teoria de 
Barenblatt, verificamos que no limite interior da zona 
onde se exercem forças de coesão, as partículas se 
encontram sensivelmente mais afastadas do que o nor- 
mal. Todavia tal afastamento não é tão grande que as 
forças de atracção se deixem de fazer sentir. Quer dizer 
que o afastamento das partículas é microscópico pelo 
menos numa parte da zona da fenda (a parte exterior) 
onde se fazem sentir as forças de coesão. 


Então o perfil da superfício da fenda apresentado 
na fig. 4 deverá ser considerado como o perfil definido 
pelas partículas situadas dum e do outro lado do plano 
central da fenda. Sendo assim a fenda apresentar-se-á, 
quando vista à escala macroscópica, com um compri- 
mento inferior ao correspondente à fig. 4, e com um 
deslocamento finito nas extremidades (macroscópicas). 
Às zonas compreendidas entre as extremidades da fenda 
vista à escala microscópica e as da fenda vista à escala 
microscópica (representada na fig. 4) encontrar-se-ão 
deformadas em parte elástica e em parte plasticamente., 


O terceiro postulado exigindo a constância (para 
um dado material em condiçes definidas) da distribuição 
das forças de coesão que correspondem à extensão da 
fanda, implica a constância do deslocamento correspon- 
dente das extremidades (macroscópicas) da fenda desda 
o início da sua extensão. Conclui-se pois que a fenda 
expandir-se-á quando aquele deslocamento à atingir de- 
terminado valor crítico,constant e para um dado mate- 


TÉCNICA 421 


rial em condições definidas, Esta conclusão concorda 
com o que já se disse a propósito do desenvolvimento 
da Wells. 

Ainda o facto, verificado experimentalmente por 
Wells, do deslocamento da abertura correspondente à 
iniciação da ruptura ser, quando esta se dá após ou 
imediatamente antes da cedência generalizada, larga- 
mente independente do comprimento da fenda, concorda 
com o terceiro postulado. 


6.5 Critério para a extensão duma fenda num 
metal! semi-dúctil 


Numa placa dum metal semi-dúctil, contendo uma 
fenda e submetida a tensões de tracção normais à fenda 
formam-se, como já se disse, nas extremidades desta 
zonas plásticas, a partir de valores relativamente baixos 
das tensões aplicadas. Ensaios efectuados em condições 
tais que as fendas se mantêm em equilíbrio estático (por- 
tanto com comprimento constante) demonstram que as 
dimensões das zonas plásticas aumentam à medida que 
aumenta a intensidade das forças aplicadas. Constata-se 
também experimentalmente que a forma das zonas plás- 
ticas é diferente consoante a fenda se encontra em equi- 
líbrio estático ou em equilíbrio dinâmico (em expansão). 
Ainda se verificou em aços que as zonas plásticas de 
fendas em equilíbrio estático têm a forma da secção 
recta da chama duma vela, 

Em certos aços e plásticos a questão tem sido re- 
duzida a termos mais simples. Assim em placas de cer- 
tos aços contendo fendas em equilibrio estático Dugdale 
constatou que as zonas plásticas são como que o pro- 
longamento da fenda (fig. 5). 


Fig. 5 


Esta constatação concorda com as constatações 
feitas sobre a teoria de Barenblatt, Dugdala formulou as 
seguintes hipóteses: 
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al Na zona plástica o material está submetido à 
tansão de cedência o 


b) àAforma e dimensões da zona plástica é tal 
que o material em regime elástico pode consi- 
rar-se limitado internamente por uma elipse 
achatada de comprimento 2 (a + p). 


c) O comprimento da zona plástica é tal que não 
exxistem singularidades de tensões nas extremi- 
dades da zona plástica. 


Por outro lado é de supor que a extensão da fenda 
comece quando a capacidade de deformação plástica 
nas suas extremidades fique esgotada. À capacidade de 
deformação plástica máxima corresponderá um valor à 
(crítico) do deslocamento da abertura da fenda. 

A partir destas premissas é possível deduzir, tal 
como fizeram Goodier e Field, que à está relacionado 


com as características do material, com o semi-compri- 
mento da fenda e com as tensões aplicadas pela seguinte 
expressão 


80 , a mo 
à mi log sec E 
q E 20, 


Desta expressão, já comprovada experimentalmeante 
em vários casos, conclui-se que quanto menor for a 
ductilidade local, de que 5 é uma medida, menor será 
a tensão o necessária à extensão da fenda, Por conse- 
quência a extensão da fenda pode dar-se sob a acção 
de baixas tensões aplicadas quando a capacidade local 
de deformação plástica for reduzida, por exemplo por 
encruamento, envelhecimento, abaixamento da tempera- 
tura, aumento da velocidade de aplicação das forças 
exteriores. 
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ANALYSIS OF FADE DURATIONS ON 
DIFFERENT MICROWAVE PATHS 


RESUMO 


Este trabalho expõe a análise feita em três feixes 
hertzianos de microondas, Essa análise conduziu à obten- 
ção de distribuições log-normais, relacionando a dura- 
ção dos fenómenos individuais de fading com o seu 
número e com o tempo total de fading. Os resultados 
são depois comparados com algumas recomendações 
do CCIR. 


1. INTRODUCTION 


This paper is concerned with the analysis of fades 
observed in propagation experiments carried out over 
three paths in Italy. With such analysis we try to deter- 
mine the log-normal distributions of single events, defi- 
ning its dependence on path characteristics (length and 
frequency). 

The paths here after analysed have already been 
studied by P, Quarta (1), [2], [3] who gave the yearly 


JOSÉ à, CC. SARAIVA MENDES 
Engenheiro Electrotécnico (IST) 


SYNOPSIS 


This paper reports the analysis of fade durations 
on three different microwave paths. We have obtained 
some log-normal distributions when relating the dura- 
tion of single events to their number and to the total 
time of fading. The results are compared with some 
CCIR recommendations, 


cumulative distributions of fading for all of them. The 
equipment used for the propagation tests allowed the 
measurement of fade durations higher than 7.5 s, with 
an error E 3 dB. 


2. PATH CHARACTERISTICS 
The three paths analysed have the following main 


characteristics: and their profiles are shown on fig. 1, 2 
and 3. 


TABLE | 
Path Type dir Er 
Beigua-Milano Over land 2000 | a 
Cimone-Milano Over land 6020 | 189 
Naples-Capri Over sea 6000 31.5 


3. DATA ANALYSIS 

All tha results shown in this paper correspond to 
measurements carried out during the summer. The ana- 
lysis corresponds to the following total durations of 


observations: 


— 923 hours for the path Beigua Milano 
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— 227 hours for the path Naples - Capri 
— 216 hours for the path Cimone - Milano 


The path Beigua - Milano (2000 MHz) has been 
observed during a longer period for it had slower varia- 
tion of the received field, as we may conclude from 
fig. 4, where we show a typical recording of the fading 
for each path. The evolution of the fading of the path 
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Naples - Capri (6200 MHz) is faster than for the path 


Cimone - Milano, (6002 MHz) for it is an over sea 
path. 


On fig. 5 we show the daily distribution of fading 
for the path Beigua-Milano. For the other two paths we 
have not made the same distribution as the influence of 
the daily period, was not so marked, The distribution 
of single events has been studied measuring their dura- 
tions at 5-10-15-20-25-30 dB bellow the value of 
free space propagation. If t is the duration of single 
events at a given level | of attenuation and p; its den- 
sity of probability, it is possible to write: 


log p: = f (pi) 
p; may be defined as follows: 


a) Number of single events with duration t, 


Total number of single events at level 1 


b) Total duration of all single events at level 1 


Total time of fading at level 1 
If: 
n:. — number of single events with duration t; 


ns — number of single events with duration t; 
um number of single events with duration t, 


We shall have accordingly to a) 


nl 
Pla = — — — 
Re tão T sas TF A, 
Na 
Pla = 3 
Ma = Ma E a =-f 
LR 
Ema = OCS IN 


BEIGUA-MILANO 


23,30 23.00 


24 


TÉCNICA 421 


in 


4.30 4.00 


CIMONE - MILANO | 


A E" lilk 
AULA TN LÊ 
rm mataram 

Lo th. 


! 


3.30 3.00 (h) 2.30 


(b) 


Fig. 5 


and accordingly to b) 


MM ti 
Pib= — 
m tt + nm beast t. 
na t; 
pp=———— 
nm t+trnma t+ + n É. 
nt 
Pmb = mm = no 


m t+ m t E na +n.t 
The distributions shown on fig. 6, 7 and 8 corres- 
pond to the integration of the density of probability 
defined in a). On the other hand fig. 9, 10 and 11 


show the distributions of the integration of he density 
of probability defined in b). Both the definitions lead to 


log-normal distributions, allthough with a higher disper- 
son for the definition b). 

These disributions are quite a good help for the 
study of a linkbut they must be referred either to a total 
number of events or to a total time of fading. 

The total time of fading can be determined by using 
tha Rayleigh type distributions. The determination of the 
total number of events is more difficult although it cor- 
responds to the best distribution. 

On the following figures, we show some laws which 
try to give a more complete description of fading, in 
order to find out some relations, making easier the use 
of the distributions of fig. 6 - 11. 

On fig. 12 we have ploted the average duration of 
fades at diferent levels, which has a logarithmic evolu- 
tion. It is interesting to notice that the plots for the 
three paths are paralell lines. 
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On the same fig. 12 we show the results proposed 
by the CCIR [4], for the average duration of fades at 
30 dB and 40 dB on a path 40 km long and at 
4000 MHz. 

Besides, the CCIR [4] proposes a law of variation 
of the average duration of fades with path characteristics 
which does not agree with the results that we have found 
for the analysed paths. (vide 5 4), 

On fig. 13 we show the analysis we have done con- 
cerning the number of events at different levels of atte- 
nuation. In the figure we plot the depth of fading against 
the number of events in percentage to the 5 dB events. 
The so obtained curve also follows a logarithmic law, 
for the considered levels of attenuation except for the 
over-sea path Naples - Capri. 

On the same figure, we also plot the results obtai- 
ned by Barsis et alt [5]. 


CET IEINL . 


TENTO 


som fl RT 

1000" are 

B00" 

600" A E 

Ke Neg o MOMIG NARA aa 

Rea itt RE HENMENE EE 
E oo IL Ni SR NE ETO MIRA Vo 
ES so AN INISUINININTI TI LA E 
é 100º Ham NUR ENNAE pe 

— EE ORE SE 

ig = RENA Hr TORO 

id PNERENS SR REITO! um. 

o O IN SN NSNEE IH 

e EI TI INN ONA SSNMINIADADARIE 

e AA NENAIS EN! Ri 


0] OM 1 2 4 E) 20 99,9 99,99 Ja 
%OF TOTAL NUMBER OF EVENTS AT EACH LEVEL 


Fig. R 


26 TÉCNICA 421 


3600" 


3000" HERRSHAAH+HHHHE Edde Tt 
aaa Das A 

voe TONE NT 

1500" PR 

Fo JE ER À, 

800” ts 

600" 

40d" 


FADE DURATIONS 
Ps 
o 
S 


FADE DURATIONS 
Rs 
o 


AU DD E E Ema Fa e Ni a. 
DO E O NIE CHE 


a ni 
or po pero 


| RES 
| EINCENCEEEEEEE 
| PRINT EEE À 


| [= Ce + Ea E E SR COSTS 
CRECI rr = AA MG o dd 
RS TRE, a é A A e a A a 


rot ENDA q E 
UI L III LAS! 


HIRISA QNMERINISa Sia 


0, 04 1 2 4 TO 20 40 50 99,9 99,99 “Jo 


%OF TOTAL TIME OF FADING AT EACH LEVTL 
Fig. 9 


NGS UNTLIN TO 


NE UNANTUIN LS 


01 04 1 2 4 10 96 99 99,9 99,99 h 


OF THE NUMBER OF EVENTS 
Fig. 10 


TÉCNICA 421 


21 


3600" 
3000" 


2000" 
1500" 


800" 


FADE DURATIONS 
ba 
S 
S 


= po ENÃ SRpEagÃa 
SS CENTER 


SAR 
TENDE 


Il ENNIO 
HERE SENSE ESA tesssta: 


ENS 
INN 
| O Ner 


NA 


OF TOTAL TIME OF FADING AT EACH LEVEL 
Fig. 11 


PERENE | mm amo inasaan aa 


FADE DURATIONS 


=p 
E E a mm mm aisais 
a E Rd E DE RR ERES 


ADM LO DOGS URLOa 


| ER 
| eee ANA do 


LE 


0,1 04 1 2 4 10 40 60 80 96 99,9 99,99 k 


28 


wOF TOTAL NUMBER OF EVENTS ME 
Fig. 12 


TÉCNICA 421 


4. INFLUENCE OF FREQUENCY, LENGTH AND TYPE 
OF PATH ON THE DISTRIBUTIONS OBTAINED 


According to the CCIR [4] the standard deviation 
of the log normal distributions of single events does 
not depend substantally on frequency and it is given 
approximately by log 2.7. Besides and actording to the 
same document the duration of fades is approximately 
inversely proportional to the length of the paths, Such 


duration is of the order of 4 to 5 seconds for 40 dB 
fades and of the order of 7 to 9 seconds for 30 dB 
fades, allways for 40 km paths at 4000 MHz. 

As the paths we have analysed are neither 40 km 
long nor at 4000 MHz, we can only use both relations 
proposed by the CCIR together, extrapolating the results 
for our paths and comparing them with experimental 
ones ,as follows: 


TABLE |l 
f (MHz) L (km Average duration at 30 dB (s) 
Path TCE DE E 
4000MHz 40 km Calculated Measured 
Cimone - Milano 1.5 4.71 1 + 1.27 217 
Beigna - Milano 0.5 3.15 4.45 => 5,7 22 
Naples - Capri 1.5 0.79 3.7 + 4.7 1.8 


We can see that the results we have obtained are 
very different of those proposed by the CCIR [4] 
except for the over-sea path Naples - Capri, However 
the results shown for the path Naples - Capri, is not a 
measured one. In fact the equipment did not allow us 
the measurements of fade durations 7.55, So the value 
at 30 dB has been caltulated assuming that the yearly 
cumulativo distribution of fading is a Raylegh-type dis- 
tribution, for attenuation levels below 20 dB. 

Considering fig. 12 it seems that the average dura- 
tion of fades depends more on the path length than on 
frequency. This result would agree with another do- 
cument of the CCIR [7] although it concerns total 
durations and not the duration of single events. 

On fig. 14 we show the influence of path length on 
the duration of fades, for the studied paths. It is interes- 
ting to notice that according to the result of fig. 14 
there would be no fading for very short paths. This has 
been confirmed through measurements carried out at 
11 GHz by P. Quarta et alt [8], [9] and [10], who 
have only found supplemental attenuation due to rain 
over two paths 10 and 20 kms long. 

From the results of fig. 13 it seems that for the 
overland paths studied in [5] and for the overland 
paths Cimone - Milano and Beigua - Milano, there is an 
influence of path length on the number of events. In fact 
for the limited number of paths analysed the number of 
avents increases with path length, 

Regarding the standard deviation of the distribu- 
tions of fig. 6-11, we cannot confirm its agreement 
with the CCIR as it reports values up to 40 dB while we 
hava measured a maximum of 30 dB. 

Nevertheless, as the distributions we hava obtained 
have a regular spacement, we may assume that the 
evolution of the distribution for 40 dB is paralell to that 
for 30 dB, with a spacement between them of the same 
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order of the spacement obtained between lower levels. 
In this case, the standard deviation here obtained would 
be similar to the value proposed by the CCIR, for the 
two overland paths, and lower than that value for the 
path Naples - Capri. 


b) Influence of the type of path 


às we have not found quoted in the literature any 
result concerning the influence of the type of the path on 
the duration of fades, and as the analysed results do not 
allow us to find acceptable relations, we only have 
verified the empirical formula: 


f 
Pa =(KX O) (q do 


Whera 


Pa — probablity of deep fades for small percentages 
of time 


K = 5.1 X 10º 
O = 0.4 (over mountain) 
= 1.0 (over plain) 
72 1 


= 

————— h = 

= Ei Ea 
F 


(over sea or coast) 


ho 02 2 
f = frequency 
d = path length (km) 


This formula has been proposed by H. Makino e 
K. Morita [11], and reported also by the CCIR [12]. 
Considering that deep fades are those in the order of 
40) dB, there is a good agreement between measured 
and calculated values. 
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5. CONCLUSIONS 


From all the measurements carried out by the 
Laboratory F. Vecchiacchi, concerning propagation tests 
we have chosen those of the paths Cimone - Milano, 
Beigua - Milano and Naples - Capri, in order to study 
tha distributions of single events. We may conclude that 
such distributions are log-normal. 

Regarding the influence of path characteristics on 
these distributions, we cannot present precise results, 
as we only have a reduced number of significant paths. 

Anyhow we have tried to carry out an analysis of 
other conclusions, comparing them with the results here 
obtained. 

This analysis have lead us to the following con- 
clusions: 


— tha relations proposed by the CCIR [4] as to 
a reduced number of paths over the United States, do 
not seem to be usable on the studied paths (vide 
table Il); 


— the standard deviation of the distributions stu- 
died over tha same paths on the United States, seems 
to have the same value that we have determined on over 
land paths; 


— it seems that other results proposed by the CCIR 
[7] for total duration of fades, are acceptable when 
applied to the duration of single events; 


— through the experimental data we have analysed 
it seems that the average duration of fades varies di- 
rectly with path length; 


— the average shape of fades when ploted versus 
time, seems to depend on the characteristics of the 
paths. So it has not been possible to define a standard 
model for its study. 
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Instalação de disjuntores ao ar livre, de 245 kV, em 
serviço na Electricité de France, com disjuntores de ar 
comprimido, tipo DLF 245, capacidade de ruptura 12 GVA, 


Fornecemos: 
Disjuntores de ar comprimido para instalações ao ar 
livre, para tensões de serviço de 72.5 até 765 kV 
Correntes nominais de 2000 até 4000 A 
Correntes de ruptura até 72 kA 

— Elevado grau de segurança 

— Construção robusta 

— Amplas possibilidades de adaptação ao 
funcionamento com potências mais elevadas de 
ruptura devido ao tipo de construção por blocos 

— Tempos minimos de disparo 

— Aptos a reengate automático rápido 


Aproximadamente 20 000 disjuntores Brown Boveri de 
ar comprimido instalados ao ar livre garantem 
em 45 países o funcionamento de rêdes de alta tensão. 


Representada em Portugal por: 
Sociedade de Electricidade 


Brown Boveri Limitada 


Rua de Sá da Bandeira, 481-2”, Porto 
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ANO XLIX 
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OBTENÇÃO DE DADOS METEOROLÓGICOS 
ATRAVÉS DE UM SATÉLITE ARTIFICIAL 


PREAMBULO PARA 3 TEXTOS 


Os 3 textos seguintes constituem respectivamente 
as 3 fases do projecto de um pequeno computador: 
programador de uma antena helicoidal, do Centro Espa- 
cial da Mulemba: 


1.º fases — «Obtenção de dados meteorológicos atra- 
vês de um satélite artificial» — Faria, 
C. B. 


2.º fase — «Tratamento, num computador IBM 711390, 
da informação ãàcerca da trajectória de um 
satélite artificialy — Levy, T.; Loureiro, A. 


3.º fase — «Projecto de um programador digital para 
para comando automático de antenas usa- 


1. PRIMEIRAS PALAVRAS 


1,1 — Um satélite artificial, vogando a uma distân- 
cia da Terra entre 200 e 36 000 km está em condições 
de ver grandes extensões do nosso globo. Nós, habitan- 
tes deste planeta, não podemos apreciar semelhante 
espectáculo, já porque estamos muito baixos, já porque 
nos encontramos mergulhados nos próprios fenómenos 
atmosféricos que pretendemos estudar. 

Inicialmente projectados para fins meteorológicos, 
os satélites americanos TIROS, ESSA, NIMBUS, ATS, 
etc., forneceram até hoje, elementos muito mais exten- 
sos, aplicados ainda à agricultura, correntes marítimas, 
formações geológicas, biologia maritima, locais arqueo- 
lógicos desaparecidos há milhares de anos, localização 
de balões-sonda, bóias marítimas, rastreio de animais 
selvagens, deslocamento de «ice-bergs», temperaturas 
da superfícis da água do mar, aperfeiçoamento de loca- 
lização de ilhas de dificil acesso, alterações da confi- 
guração das calotes polares. 

Mas é a observação directa do que se passa à 
superfície da Terra que tem merecido à ciência actual 
um maior cuidado, Não só os satélites meteorológicos 


* Eoslho Hi 
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BETTENCOURT FARIA 
Centro Espacial da Mulemba 
Luanda - Angola 


das na recepção de sinais de satélites arti- 
ficiais» — Martins, R.; Simões, |. 


Este projecto exemplar * constituiu-se como um 
dos eixos de trabalho do Grupo de Estudos de Ciberné- 
tica, da Universidade de Luanda, no ano lectivo de 
1972-1973 e, por si só, foi revelador de uma prática 
cibernética de trabalho, ausente ainda das Universi- 
dades Portuguesas. 

Embora o Centro Espacial da Mulemba não per- 
tença à Universidade de Luanda, nela desenvolve-se 
uma prática científica, que gostariamos de ver apoiada 
E neste sentido, a minha colaboração foi a mais mo- 
desta homenagem, 


Helder Coelho 


— com aperfeiçoamentos marginais —- mas outros, em 
aspecial os ERTS — (EARTH RESOURCES SATELLITES) 
que transmitem preciosas informações. 

Com a explosiva expansão, em certas áreas do 
globo, de herdades, indústrias, minas, etc. o natural 
aumento de produção de detrictos de vária espécies, 
logo se verificou a interdependência de todas as activi- 
dades humanas, & a crescente necessidade de um con- 
trole rigoroso dos mesmos. 

A Terra funciona como uma nave espacial, cuja 

auto-capacidade de manutenção é naturalmente limi- 
tada. Se se não exercer um controle no consumo, na 
natalidade, na poluição de águas e atmosfera, cedo a 
vida será impossível e o globo se transformará, como o 
seu satélite natural, a Lua, num mundo morto e sel- 
vagem, 
O controle das fontes de manutenção da humanidade 
só pode ser efectuado convenientemente desde que se 
disponha de abundante somatório de informações, sem- 
pre actualizadas, e provenientes de toda a parte do 
mundo, 


«Da organização de um curso sobre Sistemas Digitaiss — Electricidade, Maio de 1974 


Os satélites, dada a sua situação previlegiada, cen- 
tenas de quilómetros acima da Terra, são instrumentos 
da indiscutível valor para esse fim. De acordo com as 
câmaras que transportam e a distância focal de suas 
objectivas, «vêem» os diversos tipos de formações de 
terrenos ,florestas, rios, lagos, etc., isto para darmos 
apenas uma lista parcial. Se tomarmos estes mesmos 
satélites, mas equipados com sistemas que possam 
«ver» as regiões de infra-vermelhos, ultra-violetas e mi- 
cro-onda, o panorama muda radicalmente. Pode — tanta 
quanto se sabe de momento — descobrir-se se as folhas 
estão secas ou viçosas, se a terra se encontra molhada 
ou não, reconhecer-se os vários tipos de plantas, a es- 
pessura do gelo, infestação por insectos, poluição, etc, 

Muitas outras aplicações se decobriram para ustili- 
zação prática, dos instrumentos usados nos vários saté- 
lites artificiais [1]. 


1.2 — O presente trabalho, depois de uma breve des- 
crição histórica das várias tentativas e experiências 
efectuadas pelos americanos no sentido de conseguir 
satélites de observação meteorológica realmente fun- 
cionais e» seguros, é dedicado às técnicas de recepção 
de fotografias, relevo para o sector explorado no Cen- 
tro Espacial da Mulemba, em Angola, onde uma estação 
de rastreio foi projectada e construida em Julho de 1966, 
mantendo-se opsracional até à data em que este rela- 
tório é elaborado, 

Da utilidade de recepção regular destas fotografias, 
falam melhor os pedidos que àquele Centro têm sido 
feitos, tanto o Observatório de Bochum, na Alemanha 
Federal, como do TROPICAL ATLANTIC LABORATORY 
(a quem foram remetidas 300 fotografias) — e MISSÃO 
de ESTUDOS BIOCEANOLÓGICOS e de PESCAS de 
ANGOLA, que solicitou 6000 cópias de fotos, cobrindo 
uma zona compreendida entre o Golfo da Guiné e a cida- 
de do Cabo, no Atlântico-Sul. 

Atim de se fazer uma ideia da actividade média 
anual que uma estação de rastreio bem organizada pode 
conseguir, abaixo apontamos, condensamente, alguns 
alementos referentes ao «Centro» — de Luanda, 

Calcularam-se, no Centro Espacial da Mulemba, de 
Junho de 1966 a 31 de Outubro de 1971, 49,473 órbi- 
tas da satélites artificiais, assim distribuidas: 


N. de órbitras 


| 
Satélite | calculadas: 
— E = 
ESSA -2 SSI) 
ESSA —-4 3.570 
ESSA — 6 4.Iiio 
ESSA -— 8 15.960 
NIMBUS — C | 4.533 | 
NIMBUS — 3 8.950 
NIMBUS — 4 | 2.658 | 
| | 
ITOS — 1 | 5.302 | 


ATS— 3 | 490 


+ + = , —e — = 


Receberam-se, a partir destes cálculos, em cada 
ano, o seguinte quantitativo de fotografias: 


1988 — 279 | 
1967 — 591 
1968 — 1.476 | 
1969 — 1.702 
1970 2.232 
1971 — 1.559 

TEERAL. amena eins 7.830 fotos 


2. UM SATELITE ARTIFICIAL: «ESSA-8» 


2.1 Histórico 


Desde os dias, já distantes, em que o ARMY 
SIGNAL SERVICE decidiu iniciar, nos ESTADOS UNI- 
DOS DA AMÉRICA, em 1870, as primeiras observações 
meteorológicas para fins militares, até 1960 quando foi 
lançado o primeiro satélite meteorológico americano, 
o TIROS-I|, em 1 de Abril daquele ano, muito se passou 
no campo da investigação científica. Foi o início de 
uma nova era, vantajosamente beneficiada pela possi- 
bilidade de observação das camadas de nuvens a partir 
de uma altura superior a 602 km. 


Remetemos os leitores interessados na história dos 
satélites meteorológicos para o trabalho SATELITES 
ARTIFICIAIS AMERICANOS PARA FINS METEOROLÓ- 
GICOS, do mesmo autor [2], onde se encontram os ele- 
mentos detalhados sobre os mesmos. 


Apontamos aqui, como elemento de referência, O 
satélitoa ESSA-B. operacional em fins de1972, com as 
características seguintes: 


Data de lançamento: 


15 de Dezembro de 1968. 


Potência do emissor APT: 


5 watts. 


Frequência de transmissão APT: 


137.620 MHz. 


Número de linhas por imagem: 
800. 
Duração de transmissão de cada imagem: 


200 segundos. 


Intervalo de tempo entre o início de uma fotografia 
B a seguinte: 


358 segundos. 
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Fig. 1 — Aspecto típico de um satélite da série «ESSA» 


2.2 — Sistema de transmissão 


O sistema APT dos ESSA consiste numa câmara 
fotográfica, um emissor de frequência modulada, pro- 
jectado para transmitir imagens de televisão das ca- 
madas de nuvens que se encontram sob a nave espacial, 
durante as horas do dia. À transmissão é automática 
e imediata. Uma estação, em Terra, equipada com ante- 
na adequada, receptor e gravador, recebe estes sinais 
e transforma-os em fotografias, quando o satélite se 
encontra acima do horizonte, com as antenas do enge- 
nho espacial e as do posto de recepção, à vista uma 
da outra, durante o dia. 

As câmaras do satélite impressionam uma placa 
especial, qua mantém a imagem latente durante mais de 
5 minutos. À exposição dura 8 segundos que são apro- 
veitados para, no início de cada emissão, se transmitir 
um sinal de arranque, de 30 Hz e impulsos de sincro- 
nismo de amplitude correspondendo ao «branco», Estes 
impulsos automaticamente são deslocados para a es- 
querda, na fotografia, até que se fixam no extremo 
desse mesmo lado, no início de cada linha de 250 mi- 
lisegundos. No intervalo que se segue a imagem é 
pesquisada e os sinais transmitidos, formam uma ima- 
gem de 800 linhas (200 segundos), linhas estas que 
são perpendiculares ao plano de deslocamente da ór- 
bita do satélite. 


Além do sistema tradicional APT, os satélites ESSA 
transportam também equipamento DRIR (Direct Reo- 
dout Infra-Red) — leitura directa em infra-vermelhos 
— que consise num radiómetro ligado a um emissor 
transmitindo perfis em linhas de longa duração, de 
radiação emitida em reflectida, que se podem trans- 
formar em imagens de nuvans. 
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3. COMUNICAÇÃO COM O SATELITE: 
3.1 Estação receptora: 


O projecto do equipamento receptor para um emis- 
sor, de uma dada frequência e potência, suspenso no 
espaço a determinada distância, obedece a vários pará- 
metros. Se assumirmos que iremos receber sinais que 
provêm de uma fonte de energia electromagnética en- 
contrando-se, digamos, a uma altura de cerca de 10 a 
15º acima do horizonte, poderemos admitir, para uma 
distância de 1000 km, uma cobertura de cerca de 
2360 km de raio. Verifica-se, na prática, que devido a 
factores secundários, tais como refracção troposférica, 
sinais ocasionais trans-horizontes, etc., a área coberta é 
geralmente um pouco maior. 

Se, por outro lado, num mapa de projecção Mer- 
cator (meridianos e paralelos de latitude, em linhas 
rectas, cujas distâncias aumentam à medida que se 
aproximam dos polos) traçarmos as áreas cobertas por 
uma emissão proveniente do espaço, com as caracte- 
risticas atrás mencionadas, logo verificaremos que há 
deformações nas zonas polares, sendo a maior regu- 
laridade notada junto ao Equador. 

Com as fotografias são transmitidas, não sob a 
forma de imagens circulares, mas sim como quadrados, 
dá-se, portanto uma deformação trapezoida| na direc- 
ção dos polos, quando as estações se encontram em 
latitudes muito altas ou muito baixas, 

O cálculo da energia relativa, transmitida e rece- 
bida por um satélite artifical, prevê-se, a partir dos 
elementos de projecto, fornecidos pelo fabricante do 
engenho espacial, e de acordo com a existência ou não 
de um pré-amplificador de antena, como segue: 


Com sem 
pré-ampliticador pré-amplificador 


Sinal-ruído 1000 3000 1IDIDO 3000 
Distância 
Relação 33.1 dB /23.35 dB) 28.5 dB | 17.3 dB 


ESA 


Estes valores baseiam-se nos seguintes parâmetros: 


— 5 W. 
— 1,63 (2.1 dB) 


Potência emitida 
Ganho da antena do satélite 
Relação «frente-costas» da an- 


tena receptora — 20 dB (pelo menos) 


Perdas por polarização linear —3 dB 
Temperatura da antena em graus 
Kelvin — 1000º K. 
Ganho da antena receptora (re- — 12 (10.8 dB min) 
ferido a um irradiador iisotró- 
pico) — Circular (dextr.) 
Ondas estacionárias da antena 
receptora na PR AR 


Figura de ruído de pré-amplificador — 2 (3 dB] 
Centro de pré-amplificador — 23 dB 


Perdas no cabo entre a antena (com — 1 dB 
o pré-amplificador) e o receptor 

Figura de ruído do receptor 

Largura de banda do receptor 


— "0. (0 dB 
— 45 KHz 


Fazemos notar que as relações sinal-ruíido apon- 
tados acima, são consideradas como satisfatórias para 
obter fotografias com detalhe conveniente. No entanto, 
quando a distância entre a estação receptora e a antena 
é relativamente grande, digamos superior a 30 m a ins- 
talação de um pré-amplificador junto à antena é forte- 
mente aconselhável. Acrescentamos que na prática se 
verifica que a existência deste pré-amplificador melhora 
substancialmente o equipamento, porquanto fornece uma 
margem confortável para as eventuais más pontarias das 
antenas, que ocorrem quando há incorrecção de ajuste 
no equipamento de controle remoto ou cálculo pouco 
exacto, a partir de elementos de menor precisão. 


Na sua forma mais completa, uma estação recep- 
tora de sinais de satélites meteorológicos compõe-se 
do seguinte: 


a) — Antena para 136 a 138 MHz, 

bj) — Comando remoto e controle da posição da 
antena. 

c) — Pré-amplificador, 136/138 MHz. 

d) — Conversor 136/1138 — 26/28 MHz. 

e) — Receptor de frequência modulada 26/28 MHz, 

f) — Gravador de fita magnética de 2 canais. 

g) — Registador de «facsimile». 

hj — Adaptador para recepção em infra-vermelho, 

ij) — Gerador de 480 Hz, 

jj — Divisor de frequência (= 4; - 5; - 6). 

k) — Gerador de rampa, de 300 segundos, 

|) — Amplificador de audio. 

m) — Cronómetro. 

n) — Oscilocópio monitor. 

o) — Relógio. 


NOTA 1. — O registador de facsimile (alínea g) pode 
ser substituido, como veremos adiante, por 
um osciloscópio com câmara fotográfica. 


NOTA 2— O gravador de fita magnética, de 2 canais 
referido na alínea f]) —- pode evitar-se se 
não se desejam guardar as fotografias. No 
sistema osciloscópio é perfeitamente dis- 
pensável, uma vez que as fotos são feitas 
em negativo, constituindo assim material 
de fácil arquivo e referência futura, quando 
devidamente identificados. 


NOTA 3 — Nas estações que apenas desejam reteber 
imagens a preto e branco, o adaptador para 
recepção a infra-vermelho, indicado na alí- 
nea hj) — é inútil. 


NOTA 4 — Tanto o relógio como o cronómetro podem 
ter mecanismos eléctricos ou electrónicos 
(digitais). 


NOTA 5 — À estação pode funcionar sem o osciloscópio 
monitor indicado na alinea n). No entanto 
constitui um precioso auxiliar para con- 
trole constante da intensidade de imagem, 


NOTA 6 — O gerador-padrão de 480 ciclos pode ser 
substituido por qualquer outro de frequên- 
cia diferente, desde que a divisão dê, no fi- 
nal da contagem, os 4 Hz desejados para 
sincronização do equipamento, No caso 
apontado temos: 


480 : 4 = 120 
120 : 6B= 24 
24: 6B= 4 


Se dispusermos 1600 ciclos, poderiamos, 
por exemplo, fazer a divisão como segue: 


1600 : 4 = 400 
400 : 4 = 100 
100 : B= DO 
20 : B= 4 


3.2 Antena: 
Admitem-se como parâmetros desejáveis, para uma 


antena destinada à recepção de sinais de satélites arti- 
ficiais meteorológicos, os seguintes: 


Frequência — 137.000 MHz | 2 MHz 
Ganho — Suparior a 10 dB 
Largura do ângulo Plano em elevação — 50º aprox. 
recebido Plano em azimute — 70º aprox. 
Impedância — 50 ohms 
Relação de ondas 
estacionárias — 1,2 : 1 (ou melhor) 


Polarização — Circular, dextrógira 


Relação axial — 1 dB 
Relação frente-cos- 
tas - 20 dB (pelo menos) 


Existem numerosas modalidades de antenas para 
VHF, todas mais ou menos funcionais, Damos a seguir 
elementos condensados sobre duas das mais populares: 
helicoidais e do tipo Yaggi. As primeiras são preferiveis 
pela maior facilidade de montagem e menor custo, 


3.2.1 Antena Helicoidal 


É, pela relativa facilidade de construção e óptimos 
resultados, a mais indicada. 

Poda ser formada por um, dois ou mais helicóides, 
de acordo com o ganho pretendido, tendo-se em aten- 
ção que quanto maior o número de elementos tanto mais 
fechado será o ângulo de recepção. A experincia indica 
que se obtém bons sinais com um helicóide simples, 
de 6 a 8 espiras. À única desvantagem em relação a 2 
antenas, por exemplo, é que dada a sua impedância 
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— cerca de 130 ohms — necessita de um transformador 
de adaptação, de 1/4 de onda. Se usarmos cabo coa- 
xial de 50 a 75 ohms, com 2 antenas obtêm-ss directa- 
mente a impedância certa. 


Vejamos agora as características teórico-prático de 
uma antena helicoidal típica. 

Trata-se de uma hélice de modo axial que pos- 
sui um certo número de características pouco usuais. 
A mais saliente é por certo a propriedade de manter 
uma relação de radiação (e recepção) — da ordem dos 
1:1,7, com a directividade e ganho sempre na mes- 
ma direcção. 

A combinação mais conveniente para este tipo de 
antena é a que se se mostra no desenho axeno em que 
tomaremos os seguintes símbolos para descrever o 
helicoidal e seu plano de terra: 


= Diâmetro do helicoide. 

= Espaço entre espiras (centro a centro). 

= Ângulo do passo (tan-1 (S/ D)) 

= Comprimento da um espira. 

= Número de espiras. 

Comprimento axial 

Diâmetro do condutor usado no helicoide. 

Distância da primeira espira ao plano de 
terra. 

Diâmetro do plano de terra, 


| 


vorors mo 
| 


II 


É) 
| 


Se desenrolarmos uma espira e a colocarmos 
numa superfície plana, o perimetro (x D), o espaço en- 
tre espiras (8) e o ângulo do passo (a) podem expri- 
mir-se sob a forma de um triângulo. 

Condensamos a seguir o resultado que a prá- 
tica considerou favorável ao cálculo de antenas heli- 
coidais (à = comprimento de onda em metros): 


5 D = 0,32 
Ss = 0,22 
d = 0,02 (pode usar-se chapa) 
g = 0,12 
a 
- d 


A largura de banda (Lg) entre 2 pontos represen- 
tando meia-potência, é dada com grande aproximação 
por: 
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(graus) 


em que. 


À = Comprimento de onda no espaço livre. 
S = Espaço entre espiras, 

C = Perimetro de helicóide. 

n = N.º de espiras, 


|| 


O ganho em potência (G pl em relação a uma an- 
tena isotrópica polarizada circularmente, encontra-se 
usando a seguinte expressão: 


C ns 
G,=15 a 
P k À 
ou, se pretende em dB: 
C ns 
G, (dB) = 11.8 + 10 log K )º E | 


Vimos que as antenas helicoidais usadas para 
recepção de sinais de satélites artificiais tinham — 
usualmente as espiras enroladas no sentido do movimen- 
to dos ponteiros dos relógios (dextros). À vantagem des- 
ta configuração é dada pela verificação experimental de 
que a fase relativa de propagação da onda no helicoidal 
cumpre um máximo de directividade quando o condutor 
é enrolado deste modo. A relação axial (Ra) pode achar- 
-se aproximadamente por: 


2n+ 1 
in 


Ra = 


Falta-nos finalmente conhecer a impedânkia ter- 
minal (R) do helicóide que, como resultado de um 
grande número de medidas em antenas deste tipo, se 
pode avaliar, com um erro de = 15% pela seguinte 
expressão: 


Com C = perímetro do helicóide 


Se se pretende montar 2 ou 4 antenas no mesmo 
reflector (plano de terra), este terá 2,5 a 3 comprimen- 
tos de onda, de lado, e os centros dos «núclaos» das 
antenas ficarão a uma distância minima, uns dos outros, 
de 1,5 comprimento de onda, 

Os técnicos alemães, adoptam um critério baseado 
nos dados que se ilustram a seguir: 
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Analisando as medidas utilizadas por vários fabri- 
cantes na construção ide antenas helicoidais, projecta- 
das especialmente para a recepção de sinais de satéli- 
tes artificiais, entre 136 e 138 MHz (frequência central 
137 MHz), verificam-se tolerâncias que indicam não 
ser muito crítico o cálculo dos respectivos parâmetros. 
O mapa abaixo inscrito mostra o que acabamos de 
afirmar. 


| 


ROTHAMMEL | Rar s TACCO | KRAUSS 
SCHWARZ 
D 69,63 cm o 64,0 7,0 | 74 
s | 5255cm| 572 55,5 | 484 
g 28,00 cm 265 26,0 26,4 
G | 219,00 am "1600 165,0 | 176,0 


= Diâmetro do helicóide 
Espaço entre espiras 
= Distância da primeira espira ao reflector 
— Diâmetro do reflector (plano de terra). 


Doo 
| 


É de notar que eventualmente as diferenças apon- 
tadas se referem, em especial no que respeita à firma 
alomã RHODE & SCHWARZ, a uma melhoria de linea- 
ridade de resposta de frequência sem prejuizo para o 
ganho da antena, Às estruturas de suporte, — com 
excepção do helicóide, reflector, transformador e torre, 
com seus mecanismos de movimento em altura e azi- 
mute, — são totalmente feitos em material plástico, 
núcleo cenral, suportes de espiras, braçadeiras, parafu- 
sos e porcas. Nas antenas fabricadas por firmas ameri- 
canas, que tivémos ocasião de analisar em pormenor, 
só os suportes das espiras e saídas para o transforma- 
dor são isolados. 


TRANSFORMADORES DE IMPEDÂNCIA 


Vimos já como se calcula a impedância terminal 
de um helicóide. Para a frequência de 137 MHz e como 
resultado de várias medidas efectuadas, adopta-se o 


valor médio de 130 ohms. Os cabos coaxiais mais 
frequentemente usados, têm uma impedância carac- 
terística de 52 ou 75 ohms. portanto a sua ligação 
directa à antena traduzir-se-ãá numa falta de adaptação 
com importantes perdas por ondas estacionárias. 

A não ser que se instalem, num mesmo suporte, 2 
helicóides em paralelo, a montagem de um transformador 
de impedância, de 1/4 de onda, é, pois, imperativa. 

Como é geralmente sabido a relação de impedân- 
cia é dada por: 


2 = y Ze; Zs 
em que: 
Z — Impedância do transformador 
Ze — Impedância de entrada 


Zs = Impedância de saída 


Se tomarmos Zp — 13) ohms; e Z; = 75 e 55 
ohms respectivamente vem: 


Z=>=130xX75=õ98,74 ohms 
Z = 130 x 55 —84,56 ohms 


Por outro lado de impedância de um coaxial (tipo 


de transformador que escolhemos) — é-nos fornecida 
por: 

Z = 138 log 
com: 


D = Diâmetro interno do tubo exterior 
d = Diâmetro do condutor interior. 


Basta-nos, pois, escolher o diâmetro interno do 
tubo exterior (D) para se calcular facilmente o diâme- 
tro do condutor interno (id). Na prática verifica-se que 
há suficiente exactidão nas seguintes relações. 

O comprimento do transformador é de 1/4 de 
onda, o que nos dá: 


ao cabo Isolador 
coaxial 


ANTENA 


Fig. 4 
136 MHz = 2,20588 m. = À A: 4=0,55147 
137 MHz = 2,18978 m. = B B : 4 = 0,54745 
138 MHz =2,17391 m.=C C: 4=0,54348 


Adoptaremos a dimensão média, para 137 MbHz, 
arredondando, para 55 cm. 
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GANHO EM RELAÇÃO AO N.º DE ESPIRAS 


Ainda que saibamos já que quanto maior o n.º de 
espiras tanto mais elevada é o ganho da antena helicoi- 
dal, devemos ter sempre em linha de conta que dimi- 
nui ao mesmo tempo o ângulo de recepção, o que torna 
mais crítica a pontaria. Isto é importante se nos lem- 
brarmos que as próprias predições da NASA nem sem- 
pre são exactas quanto se desejaria. Vejamos, no qua- 
dro seguinte a relação verificada, na prática, entre o 
ganho, o número de espiras de um helicóide, e a aber- 
tura do ângulo de recepção (AR). 


N.º de Ganho 

espiras em dB AR 
4 Ro, 56 
5 9/0 48 
6 10,3 45 
7 11,7 42 
8 12,7 39 
9 13,0 a 
10 13,7 3E 
11 13,9 33 
12 14,0 31 


3.2.2 — Antenas YAGGI/UDA 


Quase simultaneamente no Japão e na Inglaterra, 
dois cientistas de diferente nacionalidade surgiram com 
a mesma descrição da um novo tipo de antenas; o pro- 
fessor S. Uda, da Unversidade Imperial de Tohoku 
(1926/1927) e H. Yagi (Junho 1928). Por isso mesmo 
denominada Yaggi/Uda, a antena é constituída por uma 
séria de elementos directores, dispostos linearmente, a 
que se segue o elemento excitado, atrás do qual é insta- 
lado o reflector. Tem a polarização definida pelo plano 
dos respectivos componentes, o que, na forma simples, 
a tornaria inadequada para a recepção de sinais de uma 
nave espacial cujas antenas variam de posição em relação 
à estação receptora, em Terra, 

Dada a facilidade de construção, economia de 
preço e eficiência, foi necessário estudar a sua utilização 
para a recepção de sinais de satélites artificiais, em que 
à indispensável a polarização cirbular. Obteve-se esta 
configuração mediante a montagem de duas antenas 
cruzadas ou, sendo preferivel a montagem em separado, 
num plano de 90º, inclinando uma em relação à outra, 
da 45º, 

Existe, na prática, uma variação óptima, quanto 
ao número de elementos directores e sua separação, 
afim de se obter um determinado ganho e uma melhor 
adaptação de impedância ao cabo coaxial que transpor- 
tará a energia até ao receptor (ou pré-amplificador). Por 
outro lado quanto maior o número de directores tanto 
maior também, será o ganho da antena e menor a largura 
do feixe recebido. 
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Pelo exame deste quadro se verifica que não há 
interesse em construir antenas Yaggi/Uda com mais de 
15 directores, procurando-se antes aumentar o ganho 
pela montagem de conjuntos de antenas simples ou cru- 
zadas, conforme o tipo de polarização desejado, 


Largura de 


N.º de Ganho 
directores |feixe (graus) (dB) 
3o 22 22 
20 26 21 
13 | 31 15 
9 37 13 


4 46 8 


O comprimento óptimo dos directores (em compri- 
mentos de onda), é ainda função do número de ele- 
mentos. 


Nº de Comprimento 
airasénros dos directores 
(comp. de onda) 
42 0,385 
30 0,409 
20 0,407 
13 0,414 
10 0,420 
7 0,423 
5 0,434 


Se tomarmos como mais aconselhável, na prática, 
a utilização de uma antena com 13 directores, verifi- 
camos que uma oscilação de medida destes elementos 
entra 0,5 e 10,34 de comprimento de onda, nos fornece 
uma vasta gama de impedância de entrada, adequada 
portanto ao projecto de uma antena a adaptar a deter- 
minado cabo disponível para o efeito. 

Impedância de entrada de uma antena-tipo, de 13 
directores, com reflector de um comprimento igual a 
0,5 a distância entre directores de 0,34 (medidas em 
comprimento de onda): 


Distância entre directores 
em comprimento de onda; 


Distância do reflector 
ao elemento excitado 


em comprimentos de ondas: 0,406 0,420 
0,25 62 ohms 50 ohms 
0,18 | 50 » 43 » 
0,15 32 » 27; » 
0,13 22 » mia 
; a 


0,10 | 4% 


Podemos tomar como suficientemente exactas, pa- 
ra todos os fins práticos, as seguintes relações: 
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149.000 


MHZ 
— 62.500 | 
t (MHz) 
— 142.000) 
E (MHz) 
55.000 — 
4 (MHz) 


B = 
Directores 


135.000 | 
E (MHz) 


= 132.000 
f (MHz) 


Exci tado 


m 


E. Reflector 


Pee ps 


Em que: 


Nota: medidas em rm.” 


Fig. 5 


A = Comprimento do reflector 

B = Distância do reflector ao elemento excitado 
C = Comprimento do elemento excitado 

D = Distância do elemento excitado ao primeiro 
director 

E = Distância entre directores 

F — Comprimento ido director 


Para 137 MHz, viria: 


A = 1,087 m 
B = 0,456 m 
C = 1,036 m 
D= 0,49 m 
E = 0,985 m 
F= 0,963 m 


3.3 — Pré-ameplificador 


Afim de compensar as perdas eventuais no cabo de 
ligação da antena ao receptor, é da maior conveniência 
instalar um pré-amplificador junto à antena. Trata-se 
de unidade simples que em rigor pode usar nuvistores, 
válvulas, transistores (FET's ou outros adequados para 
VHF ou UHF). Existem numerosos circuitos, todos fun- 
cionais. Modernamente preferem-se os FET's ainda que, 
nas regiões sujeitas a trovoadas não sejam muito acon- 
selháveis, pela facilidade com que uma descarga elec- 
trostática os destroi. À impedância de entrada e saida 
deve ajustar-se à da saída da antena e entrada do recep- 
tor, cuidando-se que tenha um nivel de ruído o mais 
baixo possivel, por cuidadosa neutralização. Tratando-se 
de equipamento que vai ficar exposto a intempéria, con- 
vém que esteja incluido numa caixa estanque à humi- 
dade. Alguns fabricantes utilizam invólucros de alumi- 
nio fundido com tampa fixada por vários parafusos 
e usando junta adequada. A entrada e saída de radio 
frequência, bem como os ajustes de sintonização, fa- 
zem-se por fichas instaladas em baixo, Frequentemente 
a tensão de alimentação é proveniente do próprio 
raceptor. 
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ESPECIFICAÇÕES ACONSELHÁVEIS PARA O PRÉ-AM- 
PLIFICADOR: 


Frequência central — 137 MHz 
Largura de banda a 3-dB — 5 MHz 
Figura de ruido — & dB 
Ganho — 20 a 40 dB 
Impedância de entrada — 50 ohms 


Relação de ondas estacionárias — 1.2 : 1 


3.4 — Eiltro bipolar pré-selector: 


Alguns fabricantes aconselham e adopam a insta- 
lação de um filtro pré-selector entre a antena e O con- 
versor, mas não é prática obrigatória. O filtro que me- 
lhor satisfaz o compromisso entre perdas por inserção e 
dimensões físicas, é o duplo, de cavidade, com con- 
dutor central. 


3.5 — Conversor 


O sinal recebido pela antena, seleccionada pelo fil- 
tro e subsequentemente amplificado, sofre um aumento 
de amplitude no andar frontal do conversor, antes de ser 
transformado numa frequência de 26 a 28 MHz, por 
batimento com um oscilador local, tomandado a cristal 
de quartzo. 


PARÂMETROS TÍPICOS DO CONVERSOR: 


Frequência de entrada — — 136/138 MHz 


Frequência de saida — 26/ 28 MHz 
Largura de banda — 2 MHz 
Ganho — > 2) dê 
Figura do ruído — < 2/0 dB 


Ainda que existam, modernamente, os mais varia- 
dos tipos de conversores, usando desde válvulas e nu- 
vistores, até aos semicondutores — transistores de si- 
nal, FET's, diodos tunel — etc., o princípio fundamental 
é sempre O mesmo: um dos dois andares amplificado- 
res de rádio frequência, sintonizados à frequência de 
136 a 138 MHz; o primeiro tem de ser convenientemente 
projectado e construido, pois dele depende a relação 
sinal-ruído. Há interesse em que possua alto ganho, mas 
os componentes a usar não devem contribuir, por si 
mesmos com agitação térmica demasiada. O segundo 
andar deste amplificador é acoplado ao anterior, me- 
diante um circulo duplo, sintonizado, cuja duplicidade 
estabelece a largura de banda requerida, Segue-se o mis- 
turador, também ligado ao anterior com circuito duplo, 
sintonizado. 

A saída fornece a frequência intermédia, de 26 a 
28 mHz a ser aplicada à entrada do receptor. Acoplado 
a este misturador temos o oscilador local (também por 
meio de circuito duplo) — sintonizado — que lhe in- 
jecta o sinal de mistura, e garante uniformidade de res- 
posta na frequência intermédia. Este oscilador pode 
ser auto-multiplicador ou terá um ou mais andares para 
o efeito, Em regra utiliza-se cristal de quartzo, para boa 
astabilidade, na normalidade de «overtones. 
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3.6 — Alinhamento e caracteristicas de ruido: 


O alinhamento do conversor faz-se pela seguinte 
ordem: 


a) — Oscilador local e seus eventuais andares; 

b) — Amplificador ou amplificadores de rádio fre- 
quência; 

c) — Secção de frequência intermédia, 


FIGURA DE RUÍDO: 


Usa-se um gerador de ruído de largo espectro e 
impedância igual a 5O ohms. Ligado ao conversor, a 
sua saida segue através de um atenuador comutável 
de 3 dB, para uma ponta detectora que irá indicar num 
voltimetro de corrente alterna, a saída do sinal. Compa- 
rando o nivel próprio do equipamento a plena amplitude 
— e com atenuação de 3 dB — com a que se obtém 
a partir do gerador de ruído, ficamos a conhecer a rela- 
cão que, como vimos já deve ser inferior a 2 dB. 

Os conversores utilizados no Centro Espacial da 
Mulemba, tem todos um nivel de ruído de 1,8 dB. 


3.7 — Receptor 


Naturalmente o receptor é a unidade mais com- 
plexa do sistema. Às suas caracteristicas devem, para 
o caso presente, formular-se do seguinte modo: 


Tipo — Dupla conversão, superheterodi- 
no frequncia modelada, 

Entrada — 26/28 MHz. 
Nivel de 

entrada - 90) dBm a — 50 dBm. 
Largura de 

banda — (pré-detecção) — 30 KHz 
Sensibilidade — 118 cBm ce sinal de entrada de 


rádio frequência para uma saida 
de uma relação 'sinal-ruido de 
29 dB. 


Como exemplo típico, descrevemos aseguir um dos 
3 receptores usados desde 1966 no Centro Espacial 
da Mulemba, 

Originalmente construído pela firma KAAR ENGI- 
NEERING CO, em PALO ALTO, Califórnia — USA, 
para fins de comunicações, foi modificado no C. E. M. 
para adaptação de sinais de satélites artificiais. 

Tem as seguintes características: 


Tipo FM-39X — Superheterodino, dupla conver- 
são. 
N.º de válvulas  — 15. 
Frequência de 
entrada — 26/28 MHz. 
Sensibilidade = ES NM, 


1.º oscilador — Exterior, 29/23 MHz (variável). 
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2.º oscilador — Local — 4951.666 KC. 


1H E — 5000/6000 KHz. 
2º IF — 455 KHz. 
Largura de banda — 30 KHz. 
Detecção — Discriminador. 


Saída dé audio — 8% e 600 2. 
Relação sinal-ruí- 
do -20 dB (aprox.). 


DESCRIÇÃO: 


Três transformadores de rádio frequência sintoni- 
zados a 26/28 MHz, o primeiro, simples, e os dois 
seguintes duplos, constituem o andar frontal do re- 
ceptor. 

Segue-se a 1.º misturadora que recebe o sinal am- 
plificado, proveniente do convérsor, e o mistura com o 
1.º oscilador, que no caso presente é instalado a dis- 
tância, por conveniência operacional. poder-se-ia prea- 
ver, em vez deste gerador de rádio frequência, uma 
série de cristais, oscilando na gama pretendida. Pre- 
feriu-se, no Centro Espacial de Mulemba, a primeira 
modalidade, dada não só a dificuldade de encontrar 
no mercado os cristais de quartzo para a frequência 
desejada, como ainda porque um oscilador de frequên- 
cia variável permite desde logo o ajuste manual para 
as conversões convenientes, devidas a efeitos de 
DOPPLER. 

Nos receptores do C.E.M. o cálculo de frequência 
a injectar no misturador faz-se usando o seguinte cri- 
tério: 

Vimos já que a frequência de entrada, no conver- 
sor ,era de 136/138 MHz, correspondendo a saída a 
26/28 MHz, E esta portanto a entrada no receptor. 
A primeira frequência intermédia é de 5000 a 6000 KHz, 
e a segunda de 455 KHz. Nos equipamentos adaptados 
existia já um 2.º oscilador com cristal de 4951,666 KHz 
pelo que teve de se combinar esta frequência com a do 
1.º oscilador afim de obter as frequências intermédias 
correctas. Sendo o 1.º oscilador de frequência, de tipo 
variável, era indispensável que cobrisse as frequências 
normalmente utilizadas (em APT) pelos satélites meteo- 
rológicos americanos. Até hoje foram os seguintes: 


ATE —3 — 135.600 MHz 
NOAA — 1 — 136.700 MHz 
Todos os NIMBUS — 136.950 MHz 
ESSA 2 a 6 — 137.500 MHz 
ITOS — 1 — 137.500 MHz 
ESSA -— 8 — 137.620 MHz 
| | NIMBUS | 
ATS=3 NOAA-1 | Ena 


| ESSA-0 


Fig. 6 
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O conversor (com oscilador local de 110 MHz) 
transforma estas frequncias nas seguintes: 


135.600 em 25.600 MHz 
136.700 em 26.700 MHz 
136.950 em 26.950 MHz 
137.500 em 27.500 MHz 
137.620 em 27.620 MHz 


A primeira frequência intermédia resulta da com- 
binação (soma ou diferença) deste sinal com o que for 
gerado no oscilador ajustável (1.º oscilador) e que, 
como vimos já, tem que corresponder a qualquer fre- 
quência entre os limites de 1.º FI (5000 a 6090 KHz), 
por outro lado, e como é conhecida a frequência do 
cristal do 2.º oscilador (4951,66 KHz), o batimento deste 
com a 1.º Fl terá de dar sempre a 2.º Fl, de 455 KHz. 
A frequência de injecção do OFV (Oscilador de Frequên- 
cia Variável) é pois importante, e no caso presente, em 
que foi previamente fixada a frequência do 2.º oscilador, 
calcula-se como segue: 


1º Fl=pe+ a EI 


com 


Fi — Frequênca intermédia 
Fc — Frequência do cristal 


Temos portanto, para 2.º Fl = 455 kHz 
e Fc = 4951,666 KHz: 
1.º FI = 4951,666 + 455 — 5406,666 KHz 


Subtraindo a frequência de saída do conversor, com 
a obtida nesta operação, temos finalmente a frequência 
pretendida para o oscilador variável, Para 26950 KHz 
por exemplo, seria: 


26.950.000 
— 6.406.666 


21..543.334 


mantendo-se a mesma relação de FI. 
Nas várias frequências usadas pelos satélites, temos, 


pois: 


Freg. Satélite 1.º Oscilador 


135.600 29193.334 
136.700 21293.334 
136.950 21543.334 
137.500 22093.334 
137.620 22213.334 


Donde se conclui que a gama de sintonia deste 
oscilador é muito pequena, bastando ocupar uma faixa 
entre 20 e 23 MHz, para cobrir todos os satélites APT. 

Como a maior estabilidade se consegue em fTre- 
quências mais baixas é usual utilizar o oscilador de 
comando, numa frequência 3 a 4 vezes inferior, proce- 
dendo-se à necessária multiplicação em um ou mais 
andares de rádio frequência. 


U circuito de rádiv-frequertia Fica assim cofcebido, 
do seguinte medo: (exemplo para 155.350 MHZ): 


Ê 
> 136.950 KHZ 
atire 
É] Frê-=ampliricador 


+ 136:250 Hz 


LonvErsor 
(HO MUS; 


Sa !g Ei HH Pa 
Amplificador de 
hadis-trequência 

J9HLSE kHz 1ºFT cAMists ZarT 


| Lalb. Edo 14/06+666 Lea fer gi 


L + a: E —— ds 


DO TULL? 


ElS4S SIGAHE 


2-O ans Ss 


vscilador 
fixo 


| 4951, 666kHZ 


UscIliauor 


variavel 


Fig. 7 


Seguem-se os habituais andares, limitador-discri- 
minador-audio. O discriminador neste equipamento tem 
uma largura de faixa de 50 KHz, sendo capaz de sinto- 
nizar entra 430 e 48) kHz, e está equipado com um 
micro-amperimetro dae 0-bB0 mA (zero ao centro), fun- 
damental para sintonias exactas e correcção de des- 
vios por efeitos de DOPPLER. 


3.8 — Sistema tradutor de imagem: 


Existem, em principio, dois sistemas diferentes: 
FACSIMILE e OSCILOGRAFICO. O fim de qualquer deles 
é precisamente analizar e transformar, em forma visivel 
e permanente, a informação enviada pelos satélites, 
Tanto um como o outro necessitam possuir os circuitos 
necessários para extrair a frequência de sub-portadora 
(2499 Hz) o sinal de arranque, (300 Hz) e os 4 Hz de 
sincronização de linhas. 


Especificação típica de um tradutor «Facsimile» 


Sinal de entrada — sinal de video, modulado a 2.400 Hz: 


Nivel — Q a + 20 dB 
Impedância — B0O ohms 
Sinal de saida — sinal remodulado: 
Nivel — Ajustável da0 a + 10 dB 
Impedância — 600 ohms 


Saidas de sincronização: 


Frequência — 1.2 KHz e 480 Hz 
Nivel — Ajustável de D a 4 v. 
Impedância — 500 ohms 
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Sincronismo de linhas: 


Frequência — 4d Hz 
Nível — aprox. 14 v. 
Impedância — 800 ohms 


— 35,4 linhas/cm. 
(90 linhas/polegada) 
— 240r.p.m. 
— 228,6x228,6 cm, 
(9x9 polegadas) 


Densidade de pesquisa 


Velocidade do helicóide 
Dimensão da imagem 


NOTA: Para a recepção de fotografias en infra- 
-vermelho, o equipamento deve estar preparado para 
modificar a velocidade de rotação do helicoide para 48 
linhas por minuto. 


3.9 — Equipamento oscilográftico: 


Neste sistema, as fotografias são obtidas por meio 
de uma câmara fotográfica instalada em frente a um 
tubo de raios catódicos de face plana. As linhas com 
duração de 1/4 de segundo cada uma, deslocam-se com 
regularidade, de cima para baixo, formando uma ima- 
gem completa em 220 segundos. Neste intervalo de 
tempo o obturador da objectiva mantém-se aberto, en- 
quano o diafragma foi pré-fixado para pequena abertura 
(geralmente entre f:11 ou f:16). A focagem do ponto 
electrónico é efectuada com todo o cuidado, para que 
seja do menor diâmetro possivel, e a sua imagem abso- 
lutamente circular. O controle do brilho faz-se por co- 
mando remoto afim de garantir a presença de «negros» 
nas fotografias. 

Ainda que as fotos obtidas com aparelhos de fac- 
simile, do tipo Hell, sejam um pouco melhores, quanto 
a definição, do que os conseguidos através de um osci- 
lógrafo de raios catódicos, podem os últimos conside- 
rar-se como plenamente satisfatórios para todos os fins 
práticos, tendo-se em atenção a importante economia 
resultante, quanto ao custo do equipamento. 

Por outro lado, como o processo oscilográfico utili- 
za película vulgar, (formato 6x6 — 22 a 24 DIN) não 
há problemas quanto a stock de material sensivel, es- 
pecial. 

As condições básicas, exigidas de um oscilógrafo 
de raios catódicos, para obter reproduções fotográficas 
de sinais de satélites artificias, com caracteristicas cien- 
tificamente válidas, são as seguintes: 


a) — Tubo de raios catódicos de face plana, opti- 
camente perfeita, sem graticula, FÓSFORO 
P-11, fluorescência e fosforescência azul, 


baixa persistência. 


b) — Grande estabilidade das fontes de alimenta- 
ção, a longo periodo; 


c) — Ausência de astigmatismo no ponto electró- 
nico, nas extremidades do campo; 


d) — Possibilidade de sincronismo exterior; 
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e) — Perfeita estabilidade da base de tempo à fre- 
quência de 4Hz (0,5 seg./ecm a 1 seg./ 
fem, com variação continua); 


f) — Amplificador vertical de grande linearidade, 
banda passante de corrente continua (DC) 
até cerca de 300 KHz (ImV/cm até 20 V/cm 
de sensibilidade); 


q) — Regulação de foco, intensidade, alinhamento 
a astigatismo; 


hj) ——- Cátodo do tubo de raios catódicos acessivel 
do exterior. 


3.10 — Equipamento acessório 
3.10.1 — Sincronismo 


Ainda que frequentemente o sinal de 4 Hz 
a utilizar no sincronismo do equipamento, tanto de 
facsimile como oscilográfico, seja extraído da portadora 
do satélite, no Centro Espacial da Mulemba adoptou-se 
um sistema diferente: partindo de um gerador local de 
baixa frequência — 480 Hz — faz-se a divisão até obter 
os 4 Hz. Este processo evita que, por ausência momen- 
tânea de sinal (atmosféricos, interferências, etc.) — se 
perca o sincronismo. 

A estabilidade do gerador de 480 Hz é obtida por . 
meio de diapasão, oscilando naquela frequência. À divi- 
são faz-se com multivibradores. 


e Gerador de rampa 

O deslocamento vertical do traço no écran do tubo 
de raios catódicos pode efectuar-se de diversos modos, 
admitindo que o amplificador respectivo seja de acopla- 
mento directo (D.C.) — todavia, em qualquer caso, a li- 
nearidade é fundamental porquanto qualquer perturbação 
se traduz em traços de desagradável aspecto, na reprodu- 
ção fotográfica final. No C. E. M, usa-se um gerador de 
rampa que aproveita as excelentes propriedades de in- 
tegração de um amplificador operacional, fornecendo 
não só uma rampa perfeita e isenta de anomalias, como 
ainda funcionamento cómodo, regular e seguro. 


3.11 — Marcha das operações 


Sintonizados rigorosamente os equipamentos elec- 
trónicos, afim de assegurar um máximo de sensibilidade 
a óptima relação sinal/ruido; regulado o osciloscópio — 
depois de estabilização de pelo menos, 15 minutos — 
para obter um ponto o mais pequeno possivel, circular 
e isento de astigmatismo; acertado o início e fim de 
rampa para uma imagem quadrada (num tubo de 10x8B 
em o formato deve ser, por exemplo 7x7 em — com o 
início e o fim a cerca de 5 milímetros de cada extre- 
midade); — preparadas as efemérides para pontarias 
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da antena, afim de a manter orientada para o satélite 
de minuto a minuto; as operações manuais (que em ri- 
gor podem todas ser tomadas automaticamente, sob 
comando do próprio satélite, que, como vimos, trans- 
mite, no início de cada imagem, um sinal de comando, 
de 300 Hz) — são as seguintes: 


a) — Verificação da intensidade do traço no tubo ' 


de raios catódicos, de modo a que, nem seja 
tão intenso que «queime» a película, — ori- 
ginando o desaparecimento de pormenor — 
nem tão fraco que obrigue a não ficarem 
registados os «negros» da fotografia. 


b) — Logo que seja recebido o sinal de arranque, 
dar início à rampa linear. 


c) — Deslocar o sinal de sincronização de modo 
a que se fixe na extremidade esquerda ou 
direita da imagem. 


d) — Abrir o obturador da câmara fotográfica que 
se manterá por si mesmo aberto, na posição 
de «póser» — enquanto, linha-a-linha a foto- 
grafia se for formando, 


e) — Termonada a transmissão da informação foto- 
gráfica — reconhecida pelo aparecimento de 
um sinal contínuo de 24009 Hz — fechar o 
obturador e rodar o maquinismo, que porá 
nova porção de película sensivel em frente 
da objectiva. 
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ANO XLIX 


VOLUME XXXVI 


ENERGIE (*) 


TECHNIQUES DE L'AVENIR 


La situation de | ênergie dans 10 ou 20 ans dépen- 
dra das décisions prises aujourd'hui en ci qui concerne 
da recherche et le dévevloppment, ainsi que la construc- 
tion de centrales énergétiques. 


Ce fait a servi de base débats d'une récente réunion 
qui a eu lieu à Genêve et durant laquelle des délégués 
de 14 pays ont étudié I'évolution de nouvelles techniques 
de l'énergie, larriêre-plan économique et écologique 
aux innovations dans le secteur de I'éênergie, et le róle 
de la Comission pour |'Europe de | Organisation des Na- 
tions Unies (CEE) ia promotion de l'innovation en ma- 
tiêgre d' énergie, 


Cette réunica a eu ligu sous les auspices des Von- 
seillers des gouvernements des pays de la CEE pour 
la spience et la technique, dont une des activités prin- 
cipales est dédiée aux innovations dans le domaine de 
Vénergie. 


Soulignant la durée de chaque étape de I'évolution 
des techniques énergétiques, une étude préparée pour 
la réunion par le Secrétariat de la CEE mentionne | 'exem- 
ple du réacteur surrégénéraeur à neutrons rapides (FBR). 
Les recherches sur le FBR ont débuté vers la fin des 
années 40. Plusieurs petits prototypes expérimentaux ont 
été construis vers la fin des années 50D — par exemple 
à Dounreay en Ecosse e en URSS. Vers la fin des 
années 196) et au début des années 1970, des proto- 
types plus conséquents ont été construits afin de démon- 
trer les possibilitês commerciales du FBR. Si l'on compte 
deux ans pour les essais sur de grandes installations de 
démonstration, et cing à sept ans afin de concevoir et 
construire les FBR commerciaux, ce genre de centrale 
nucleáire pourrait être commercialisé aux environs de 
1985 — soit quelque 35 ans aprês le début des recher- 
ches. 


COUP D'OEIL SUR UN AVENIR PROCHE 


Les limites restreintes dans lesquelles I|'électricité 
peut immédiatement remplacer le combustible liquide en 
tant que source d'énergie dans le le secteur des trans- 
ports montre bien que du temps sera nécessaire pour 
effectuer le changement d'une technique à l'autre, || 
existe toutefois plusieurs techniques offrant d'intéressan 
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tes possibilitás à court ou à moyen terme. Quelques- 
-unes d'entre elles ont vu le jour il y a longtemps mais, 
pour des raisons économiques ou autres, nont jamais 
été largement utilisées, 


Une technique intéressante pour s attaquer aux 
problêmes de l'offre d'offre d énergie à court terme 
consiste en la gazéification et la liquéfaction du charbon; 
la gazéification semble être un moyen économique pour 
|'élimination des éêmissions de soufre qui crée un sérieux 
handicap vis-à-vis de Ievironnement. 


Des recherches sur la gazéification se poursuivent 
dans plusieurs pays afin d'obtenir un gaz à pouvoir 
calorifique élevé, sans pour cela subir les inconvénients 
de la méthode allemande Lurgi, ancienne déjã d une 
quarantaine d années. On espére que la commercialisa- 
tion pourra entrer en vigueur d'ici deux ans, 


D importantes recherches sur la gazéification sur 
place ont été effectuées en Pologne et en URSS depuis 
la début des annés 1960. 


La conversion du charbon en une huile combustible 
a atteint le stade des installations pilotes aux Etas-Unis. 
Un des problêmes techniques posés par le charbon li- 
quéfié concerne l'éliminaion de l'oxyde d'azote qu'il 
contient, afin d assurer la compatibilité avec les normes 
de pollution. 


Entretemps, le charbon entre pour une large part 
dans les plans de construction des centrales thermiques. 
En URSS par example, la production d'ênergie électrique 
sera augmentée dans l'avenir immédiat principalement 
la construction de grandes centrales sur l'emplacement 
de riches gisements charbonniers dans les bassins d'Eki- 
bastuz et dans certains autres secteurs. Les Etats-Unis 
et certains pays d'Europe occidentale, qui avaient ralen- 
ti ou diminué la production de leur industries miniêre, 
envisagent aujurd'hui da développer lutilization du 
charbon. 


SCHISTES BITUMINEUX ET SABLES ASPHALTIQUES 


Dans certains pays tels que le Canada et les Etats- 
“Units qui prossédent de larges gisements de schistes 
bitumineux et de sables asphaltiques, la commercialisa- 
tion pourrait être réalisable d'ici le début des années 


1980 à la solution des problêmes techniques et écologi- 
ques associés à la conversion de ces matiéres primigres 
en combustible. 


EVOLUTION NUCLÉAIRE: AVANTAGES ET INCONVE- 
NIENTS 


Les centrales nucléaires de la «premigre génération» 
ces avantages économiques sur les centrales thermiques 
conventionnelles, avantages qui augmentent en même 
temps que les prix des combustibles fossiles. 


Les réacteurs nucléaires actuels utilisent moins de 
un pour cent de | énergie naturelle de I'uranium et leur 
rendement thermodynamique est seulement d'environ 
32%. Les réacteurs refroidis au gaz pourront améliorer 
ce rendement, mais, si Ion on découvvre pas de plus 
grandes quantités d'uranium, une efficacité réduite cons- 
tituera un désavantage sérieux pour les centrales nuclé- 
aires dans moins de 15 ans. 


En |" absence de telles découvertes ou d'améliora- 
tions de la technique, le monde pourrait avoir à faire 
face à la même pénurie d'uranium avant l'an 2000 que 
celle prévue à l|'heure actuelle pour les combustibles 
fossiles. 


Toutefois, la course à I'énergie de «premiêre géné- 
ration» continue, On estime que, d'ici 1980, 25 pour 
cent de |'énergie életrique et 10 pour cent de I'énergie 
totale produits aux Etats-Unis proviendront de centrales 
nucléaires. Dee 1974 à 1987, on estime que | Europe 
occidentale consacrera environ 120 milliards de dollars 
à la construction de centrales nucléaires. LURSS exé- 
cute actuellement un programe de construction de cen- 
trales nucléaires, échelonné sur 10 ou 12 ans, et des 
réacteurs de la «premiére génération» continuent à être 
mis au point. Des centrales nucléaires fonctionent en 
Tchécoslovaquie et en République démocratique alle- 
mande, et d'autres vont vont bientôt être mises en ser- 
vice en Bulgarie, en Hongrie, en Pologne et en Rouma- 
nie. Cette tendance va se poursuivre jusqu à ce que les 
réacteurs surrégénérateurs à fission et les techniques de 
fusion soint commercialisés. 


ENERGIE GÉOTHERMIQUE 


Bon marché, et de production relativement simple, 
I'énergie géothermique ne joue toutefois qu'un rôle mi- 
neur sur la scêne énergétique mondiale, bien que la 
premiêre centrale géothermal ait été installée il y a 70 
ans à bLardarelllo en lItalie. La Nouvelle-Zélande, les 
Etats-Unis, le Japon et I'URSS — les plus importants 
exploitants d énergie géothermique — ont suivi I'talie; 
la Mexique est en train de construire une centrale. Il 
existe également des réservoirs moins chauds qui sont, 
ou peuvent être, utilisés en remplacement de lI'énergie. 
La capitale de l'Islande, Reykjavik, est presque entiêre- 
ment chauffée à l'eau chaude d'origine géothermique, 
ainsi que certaires villes de Hongrie ea d' URSS. 


Les progrês techniques permettent une exploitation 
plus poussée des ressources géothermiques. Autrefois 
d'un usage três localisé, lénergie géothermique peut 
maintenant approvisionner des centrales éloignées grãce 
aux techniques de transmission par haut voltage. L'hydro- 
gêne sulfureux, qui constitui un problême ecologique lié 
aux opérations géothermiques, peut être extrait chimi- 
quement du gaz non condensable e converti en soufre 
élémentaire, Enfin, de nouvelles techniques de forage 
et de complétion des puits pourraient être la source de 
possibilités de développement plus étendues en ce qui 
concerne lI'énergie géothermique. 


SOURCE D'ENERGIE À PLUS LONG TERME 


Entre autres sources d'énergie à long terme, trois 
— les réateurs surrégénérateurs (FBRs), les réacteurs 
de fusion nucléaire et I'énergie solaire — présentent des 
perspectives prometteuses. 


Les FBRs constituent le point d'aboutissement de 
l'évolution des techniques de fission. lis permettront de 
récupérer 60 à 75% de toute |'énergie disponible émise 
par Iuranium utilisé comme combustible et de créer plus 
de matériaux fissiles qu'ils n'en utilisent. On prévoit que 
les FBRs ne consommeront que 1,3 tonne duranium 
par million de kW/an d électricité produite, contre 171 
tonnes pour les réacteurs de «prmiêre génération», En 
1973, VURSS a mis en exploitation avec sucês un pro- 
totype industriel du FBR. Un second est actuellement en 
construction, destiné à la centrale de Biéloyarsk, dans 
[Qural, Au début de 1974, le Phénix français a été le 
premier FBR d'Europe occidentale à être porté à sa 
puissance maximale. Le projet trinational Super-Phénix, 
comprenant la France, la République fédéral d'Allema- 
gne et I'Italie, pourrait commencer à fonctioner au début 
des années 1980, étant entendu que 5 ans seront né- 
cessaires à sa construction, La République fédérale 
d'Allemagne, le Japon, le Royaume-Uni et les Etats-Unis 
construisent des FBRs de démonstration, 


Deux techniques de contrôle de la fusion sont 
actuellement à l'étude. Elles suscitent beaucoup d'es- 
poirs en vue de la solution du problêms des besoins 
d'énergie mondiaux à long terme, L'étude appuvye l'idée 
que des recherches parallêles sur deux procédés diffé- 
rents visant au même but pourraient présenter des 
avantages en termes de coút et de temps; elles consti 
ent une assurance contre |'échec, et lorsque les deux sol- 
utions s'avêrent réalisables, elles permettent de choisir 
la plus économique d'entre elles. 


Une de ces techniques vise à obtenir la fusion 
thermonucléaire en élevant la température de |I'hydro- 
gêne servant de combustible 3 environ 1 million de 
degrés centigrades à l'intérieur d'un «conteneur» magné- 
tique. La seconde consiste à mitrailler une boulette de 
deutérium et de tritium par de nombreuses décharges 
de lasers qui la portent à des températeures assez éle- 
vées pour amorcer la fusion. Ce processus, qui élimine 
la problême du confinent magnêtique, utilise une en- 
ceinte sphérique remplice de lithium fondu qui absorde 
la chaleur. Lénergie puissante obtenue par une série 
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continue de micro-explosions thermonucléaires est trans- 
mise au Tithium fondu environnant qui, à son tour, est 
utilisé afin de produire la vapeur alimentant une centrale 
électrique. 


Dés recherches sur ces deux techniques sont en 
cours aux Etats-Unis et en URSS, dans |'espoir que, 
vers lan 2000, les réacteurs à fusion nucléaire pourront 
être commercialisés, 


Les réacteurs de fusion possêdent de nombraux 
avantages sur les réacteurs de fission, entre autres la 
production d'une plus grande quantité d énergie pour 
chaque unité de combustible, En outre, il est vraisem- 
blabe que le coút d'investissement pour une centrale de 
fusion sera moins élevé que pour un réacteur surrégé- 
nérateur. Les réacteurs de fusion causeront moins de 
pollution, car ils n émettent pas de caesium ou de stron- 
tium radioactifs, et il est impossible qu'une centrale 
cevienne incontrôlable et fonde ou explose. Si le procêdé 
de fusion se déréglait, la réaction stopperait instanta- 
nément. 


Les systêmes utilisant l'énergie solaire pour le 
chauffage et la climatisation des logements ont donné 
de bons résultats au stade de la démonstration et 
d'autres aplications pourront être commercialisées dans 
un délai de cing ans. Daprês l|'êtude, ces systêmes 
pourraient devenir un moyen important de conservation 
de lapprovisionnement en énergie. Mais elle souligne 
que les investissements de bas sont élevés et que si un 
moyen pratique d'emmagasiner la chaleur solaire nest 
pas trouvé, il faudra recourir à d'autres sources 
dénergie afin de soutenir le chauffage solaire 
pendant la nuit et les longues périodes de mauvais 
temps. L'embloi à grande échelle de lIénergie solaire 
pour la production d électricité ou de combustibles se- 
condaires risque donc d' appartenir encore à un avenir 
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assez éloignê. Une des possibilitês consisterait à accu- 
muler | énergis solaire dans des stations sattelites et de 
la transmetre à des sations au sol, Une autre implique 
la création de zones collectrices d énerggie solaire au 
sol comme moyen d'utiliser la lumiêre directe du soleil 
pour produire de la vapeur, à son tour génératrice d éner- 
gie électrique, 

Le charbon continuera de figurer, même à long ter- 
me, parmi les principales sources d énergie; la magné- 
tohydrodynamique (MHD) est une technique nouvelle et 
prometteuse d utilisation efficace et propre du charbon. 
Elle offre la possibilitê de convertir la chaleur. des gas 
de combustion directement en electricitê. Les centrales 
MHD pourraient atteindre un rendement de 50, ou mê- 
me 6) pour cent, comparées aux centrales actuelles à 
combustibles fossiles (40%). Elles pourraient égale- 
ment constituer l'une des meilleures méthodes d'élimi- 
nation des oxydes de soufre, et de réduction des émis- 
sions d oxydes d'azote des centrales chaufées au char- 
bon; elles produiraient également moins de pollution 
thermique. 


Les cenirales nucléaires et solaires de lavenir 
produiront de lI'énergie électrique, mais elles pourront 
également utiliser leur énergie pour créer un gaz dhy 
drogêne — combustible synthétique secondaire présen- 
tant de nombreux avantages économiques et écologiques, 
L'hydrogêne peut être produit aisément et emmagasiné 
économiquement. Lorsquil a brúlé, son seul résidu est 
l'eau. Le transport demeure le principal problême. Dans 
une future «économie de Ihydrogêne», les pipelines 
religraient les grandes centrales nucléaires ou solaires 
aux ligux de consommation de ce gas oú ce dernier 
serait brúlé afin de produir de I'électricité, de la chaleur, 
ou de | énergise mécanique. 
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Synopsis of articles published in «Técnica» n.º 42] 
XLIX — October 1974 


U. D. C. 621,392.1:621,396.029.1 
JOSE 4, C. SARAIVA MENDES 


ANALYSIS OF FADE DURATIONS ON DIFFERENT 
MICROWAVE PATHS 
Técnica N.º 421 — XLIX — 10-1974, p. 23 - 32 
This paper reports the analysis of fade durations on three 
different miorowave paths. We have obtained some log- 
-normal distributions when relating the duration of fading. 


The results are compared wlh some CCIR recommen- 
dations. 


U. D. €C. 551.501:629.19 
BETTENCOURT FARIA 
METECROLOGIC SATELLITES 
Técnica N.º 421 — XLIX — 10-1974, p. 33 - 44 


U. D. C. 624.04:518.6 


E. R. DE ARANTES E OLIVEIRA 
L. F. GONÇALVES DE AZEVEDO 


THE CONJUGATED GRADIANT METHOD AND ITS 
APLICATION IN STRUCTURAL ANALYSIS 
Técnica 421 XLIX — 10-1974, p. 1-6 

lt is shown that the use of the conjugated gradiant techni- 

que in the solution of the systems of equations associated 

with tha analysis of structures makes it possible to work 
directly wits the construction of the global stifness matrix. 


U. D. C. 539.31:518.6 


E. R. DE ARANTES E OLIVEIRA 
R. TEIXEIRA DUARTE 


à FICTITIOUS FORCE METHOD FOR THE SOLUTION 
OF GEOMETRICALLY NON-LINEAR PROBLEMS 
Técnica Nº 421 — XLIX— 104974, p. 7-12 

à method for the solution of geometrically non-linear pro- 
blems is presented which is characterized by the use of 
fictitious forces obtained without resorting to integrations 


over tha elements or to the non-linear form of he strain- 
“displacement equations. 


SALGADO PRATA 
INTRODUCTION TO THE MECHANICAL FRACTURE 
Técnica 421 — ALIX — 10-1974, p. 13-21 


U. D. C. 539.55 


The fracture mechanic relates the properties, the compacity 
faults of materials and its condition of aplication with 
the fracture phenomenon. 

This work presents basic concepts and analytic for- 
mujations of the fracture mechanies. 
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Resumo dos artigos publicados na «Técnica» n.º 42] 
e Ano XLIX — Outubro 1974 


Cc. D. U, 624.04:518.6 


E. R. DE ARANTES E OLIVEIRA 
L. F. GONÇALVES DE AZEVEDO 


O METODO DOS GRADIANTES CONJUGADOS 
E à SUA APLICAÇÃO NA ANALISE DE ESTRUTURAS 


Tóznica N.º 421 — XLIX — 10-1974, p. 1-6 


Mostra-se que a utilização do método dos gradiantes con- 
jugados na resolução dos sis'emas de equações associadas 
com a análise de estruturas torna nossivel trabalhar directa- 
menta com a matrizes de rigidez dos elementos e, deste 
modo, evitar a formação da matriz de rigidez global. 


C. D. U. 539.31:518.6 


E. R. ARANTES E OLIVEIRA 
R. TEIXEIRA DUARTE 


UM METODO DE FORÇAS FICTICIAS 
PARA A RESOLUÇÃO DE PROBLEMAS 
GEOMETRICAMENTE NÃO LINEARES 


Tócnica N.º 421 — XLIX — 10-1974, p. 7-12 


Apresenta-se um método para a resolução de problemas 
geometricamenta não lineares caracterizado pela uiilização 
de forças fictícias obtidas sem qualquer recurso a integra- 
ções sobre os elementos ou à forma não-linear das equa- 
cões deformações-deslocamentos. 


SALGADO PRATA C. D. U. 539.66 
INTRODUÇÃO A MECANICA DA FRACTURA 
Técnica N.º 421 — XLIX — 101974, p. 13-21 


a mecânica da fractura relaciona as propriedades, defeitos 
de compacidade dos materiais e suas condições de apli- 
cação com o fenómeno da fractura. 

Nesta trabalho apresenta-se conceitos e formulações ana- 
líticas básicos da mecânica da fractura. 
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ANALISE DA DURAÇÃO DO AMORTECIMENTO 
EM DIFERENTES FEIXES DE MICROONDAS 


Técnica N.º 421 — XLIX — 10-1974, p. 23-32 


Este trabalho expõe a análise feita em três feixes hertzianos 
do microondas. Essa aná[ise conduziu à obtenção de dis- 
tribuições log-normais, relacionando a duração dos fenó- 
menos individuais de fading com o seu número e com o 
tempo total de fading. Os resultados são depois comparados 
com algumas recomendações do CCIR. 
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